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ABSTRACT 


The  integrity  of  brittle  elements  subjected  to  time-dependent 
thermal  and  mechanical  loads  is  predicted  by  means  of  an  algorithm 
involving  three  steps:  the  determination  of  the  temperature  dis¬ 
tribution,  the  determination  of  the  thermal  and  mechanical  stresses, 
and  the  statistical  description  of  the  resistance  of  the  brittle 
material.  Experiments  are  conducted  on  beams  and  disks  under  tran¬ 
sient  thermal  and  mechanical  loads.  The  fracture  probability-time 
curves  obtained  for  these  members  are  very  accurately  predicted  by 
the  fracture  algorithm.  The  sensitivity  of  the  beam  and  disk  re¬ 
sponse  to  changes  in  the  elastic,  thermal,  and  statistical  strength 
parameters  is  investigated.  A  combined  stress  fracture  theory  for 
brittle  elements  is  developed  which  accounts  for  a  history  of  ther¬ 
mal  and  mechanical  loading. 
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SECTION  I 


INTRODUCTION 

Culminating  several  years  of  research  in  the  application  and 
development  of  statistical  fracture  theory,  this  report  addresses 
itself  to  the  problem  of  predicting  the  structural  integrity  of 
ceramic  elements  which  are  subjected  to  transient  thermal  and  me¬ 
chanical  loading.  Our  specific  objectives  were  fourfold: 

(1)  Develop  an  analysis  procedure  for  the  "thermal  shock 
response"  of  brictle  materials. 

(2)  Verify  the  analysis  procedure  using  simple  ceramic 
elements . 

(3)  Develop  a  combined  stress  theory  for  brittle 
materials. 

(4)  Describe  the  inadequacies  or  shortcomings  of  the 
analysis  procedure. 

For  reasons  that  will  be  made  clear,  we  were  unable  to  respond  to 
the  latter  objective. 

Before  we  describe  our  general  method  of  attack,  a  few  corn- 
merits  appear  to  be  in  order  concerning  the  term  "thermal  shock." 
The  sudden  cooling,  which  takes  place  when  a  very  hot  solid  is 
plunged  Into  a  relatively  cool  liquid,  creates  tensile  stresses 
on  the  surface  of  the  solid  which  frequently  cause  cracking.  The 
expression  "thermal  shock,"  which  is  usually  associated  with 
this  rapid  cooling  process,  seems  to  he  sort  of  a  biological  des¬ 
cription  of  the  solid's  reaction.  In  the  jargon  of  the  physical 
sciences,  the  term  shock  has  come  to  represent  a  stationary  or 
propagating  discontinuity  across  some  surface  in  a  continuum. 

This  is  not  what  is  experienced  by  a  quenched  solid  which  does 
not,  for  example,  develop  a  shock  wave.  The  complications  which 
attend  the  study  of  shock  conditions  fortunately  do  not  arise  in 
the  rapid  heating  or  cooling  of  solids.  Furthermore,  the  heating 
rates  normally  associated  with  leading  edges  and  nose  cones  are 
of  a  lower  order  of  magnitude  than  those  associated  with  quench¬ 
ing. 
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The  behavior  of  brittle  bodies  under  transient  thermal  load¬ 
ings  has  been  treated  historically  as  a  distinct  material  property 
called  thermal  shock  resistance.  With  continuing  study  it  has 
become  increasingly  clear  that  such  behavior  is  really  a  composite 
of  more  basic  physical  phenomena.  In  particular,  the  broad  prob¬ 
lem  of  strength  prediction  under  thermal  loading  can  almost  be 
viewed  as  a  classic  problem  in  continuum  strength  analysis.  It 
is  now  customary  to  divide  the  problem  into  the  following  parts: 

(a)  Determination  of  the  thermal  and  mechanical 
boundary  conditions 

(b)  Determination  of  the  time -dependent  temperature 
distribution 

(c)  Determination  of  the  time- dependent  thermal 
stress  distribution 

(d)  Assessment  of  the  temperature-dependent 
strength  or  resistance  of  the  body. 

In  the  first  three  cases  the  physical  description  of  the 
processes  involved  and  the  methods  of  approach  are  well  known; 
however,  there  are  many  analytical  and  computational  difficulties 
which  make  the  general  problem  very  -complicated.  On  the  other 
hand,  the  latter  case  is  concerned  with  a  failure  theory  for  brit¬ 
tle  materials  and  the  physics  of  this  problem  becomes  a  research 
subject,  quite  apart  from  any  possible  computational  difficulties. 
For  this  reason,  we  have  directed  most  of  our  attention  to  the 
strength  aspect  of  the  response  problem. 

An  experimental  program  was  designed  to  minimize  the  uncer¬ 
tainties  which  might  arise  in  the  temperature  and  stress  determ¬ 
inations.  For  example,  rather  than  calculate  the  temperature  dis¬ 
tribution  and  expose  ourselves  to  possible  imprecision  in  the 
characterization  of  the  boundary  conditions,  we  measured  the  dis¬ 
tribution  with  thermocouples.  To  reduce  the  thermal  and  mechanical 
stress  analysis  to  the  simpLest  possible  form,  elements  were 
sought  for  which  a  one-dimensional  stress  analysis  was  possible; 
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we  choose  a  beam  subjected  to  terminal  couples  and  heated  along 
the  bottom  surface  and  a  circular  disk  which  was  heated  around 
the  inside  edge  of  a  concentric  hole. 

A  high  purity  alumina  was  selected  for  our  experiments  be¬ 
cause  of  the  wide  experience  accumulated  with  such  materials.  This 
experience  manifests  itself  in  a  consistent  manufacturing  capabi¬ 
lity  and  in  plentiful  data  for  both  mechanical  and  thermal  prop¬ 
erties.  Some  60  beams  and  60  disks  were  fabricated  using  Wesgo 
A1  995  and  each  member  was  exposed  to  a  time  varying  temperature 
input.  The  beam  was  subjected,  in  addition,  to  monotonically 
nondecreasing  terminal  couples. 

The  time  to  failure  for  every  test  was  recorded  and  the  dis¬ 
tribution  of  these  times  was  displayed  for  each  element  by  a 
cumulative  distribution  curve.  The  resulting  two  curves  described 
the  tradeoff  between  the  fracture  probability  and  the  failure 
time  for  the  beam  and  disk.  The  analytical  prediction  of  these 
curves  required  as  input  the  mechanical  loading  and  the  tempera¬ 
ture  distribution  as  functions  of  time.  The  following  scalar 
functions  of  temperature  were  also  required:  modulus  of  elasticity, 
Poisson's  ratio,  the  thermal  strain,  and  the  statistical  strength 
parameters  for  pure  tension.  These  curves  were  integrated  into 
a  statistical  thermal  and  mechanical  response  algorithm  which 
embraced  the  conservative  assumption  that  alumina  is  a  series  or 
weakest  link  material.  It  was  further  assumed  that  the  tensile 
distribution  curves  for  the  temperatures  of  interest  do  not  depend 
explicitly  on  time  or  load  history,  so  that,  the  effects  of 
creep,  static  fatigue,  or  general  fatigue  are  not  taken  into 
account. 

The  original  program  strategy  required  that  all  of  the 
temperature  dependent  properties  be  obtained  from  the  literature 
with  the  exception  of  the  statistical  strength  parameters  which 
we  wanted  to  develop  ourselves.  Unfortunately,  the  tension  mem¬ 
bers  ordered  for  this  purpose  were  badly  warped  and  had  to  be 
disregarded.  Since  it  was  not  expedient  to  correct  or  reorder 
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the  tension  specimens,  we  set  out  to  find  the  "tension  behavior" 
that  would  predict  the  fracture- time  distributions  that  were  ob¬ 
tained  experimentally. 

Indeed,  we  found  a  set  of  Weibull  parameters  which  made  our 
predictions  for  both  the  beam  and  the  disk  coincide  almost  ex¬ 
actly  with  the  measured  results.  This  seemed  remarkable  con¬ 
sidering  that  the  beam  stresses  were  uniaxial  and  mostly  mechanical 
and  that  the  disk  stresses  were  biaxial  and  entirely  thermal. 
Furthermore,  the  tensile  strength  distribution  that  we  assumed 
was  similar  to  the  stronger  results  described  in  the  literature. 

While  this  report  was  being  prepared,  we  obtained  a  set  of 
tension  data  from  Southern  Research  Institute  that  was  generated 
for  a  99.5  percent  alumina  using  a  gas  bearing  machine.  Our  as¬ 
sumed  tension  distribution  fits  this  data  precisely. 

The  agreement  between  the  theoretical  and  experimental 
results  for  the  beam,  disk  and  tension  specimen  is  so  close  that 
it  precludes  the  error  analysis  alluded  to  in  our  fourth  objective. 
Instead,  we  have  presented  a  variation  of  parameters  study  which 
illustrates  the  effects  on  the  fracture- time  curve  of  a  +  10  per¬ 
cent  variation  in  any  physical  property  used  in  our  prediction 
scheme. 

In  view  of  the  remarkable  results  obtained  in  this  investi¬ 
gation,  we  feel  justified  in  claiming  to  have  established  the 
potential  of  the  statistical  fracture  theory  approach  to  the 
thermal  shock  problem.  It  should  be  clearly  understood  that  es¬ 
tablishing  a  methodh  potential  and  "proving"  its  applicability  are 
very  different  accomplishments . 


SECTION  II ‘ 

THERMALLY  AND  MECHANICALLY  INDUCED  FRACTURE 


1.  INTRODUCTION 

Procedures  for  predicting  the  performance  of  a  structure 
composed  of  a  brittle  material  have  been  presented  in  our  pre¬ 
vious  work  (References  1  and  2) .  In  this  section  these  proce¬ 
dures  are  extended  to  account  for  time  varying  mechanical  and 
thermal  loading  and  temperature  dependent  material  properties. 
The  general  concept  of  a  generalized  stress  is  replaced  by  an 
approximate  biaxial  Weibull  type  formulation  for  the  "risk  of 
rupture."  This  is  a  more  restrictive  statement  and  that  this 
is  a  proper  assumption  must  be  justified  in  each  application 
of  the  procedure.  Also,  new  attention  is  drawn  to  the  role  of 
the  volume  of  the  basic  tensile  strength  distribution  specimen 
and  to  th^  fact  that  care  is  necessary  to  ensure  that  the  pro¬ 
cedure  is  always  yielding  at  least  conservative  predictions. 

2.  ASSUMPTIONS 

The  applicability  of  the  statistical  analysis  procedure 
or  fracture  algorithm  in  a  given  situation  rests  entirely  upon 
the  accuracy  of  the  basic  assumptions  employed  in  the  algorithm. 
For  this  x'eason  it  is  desirable  to  state  and  examine  all  the 
assumptions  that  enter  into  the  fracture  algorithm. 

It  is  assumed  that  the  probability  of  failure  of  a  sub¬ 
volume  of  a  brittle  structure  depends  only  upon  its  temperature, 
state  of  stress  and  volume.  Effects  such  as  creep,  strain  rate 
dependance,  stress  gradient  dependance,  static  fatigue  and 
cyclic  fatigue  are  assumed  to  be  insignificant  if  not  entirely 
absent.  Surface  effects  such  as  surface  finish  are  assumed 
to  be  insignificant  and  thus  the  strength  of  a  subvolume  does 
not  depend  upon  whether  or  not  it  is  located  on  the  surface  of 
the  structure.  All  the  material  in  the  basic  tension  specimens 
and  in  the  brittle  structure  being  analyzed  is  assumed  to  come 
from  the  same  statistical  population.  The  mechanical  and 
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thermal  loadings  are  assumed  to  be  known  deterministically. 

Also,  in  the  stress  analysis  of  the  brittle  structure,  it  is 
assumed  that  the  material  is  homogeneous  at  uniform  temperature, 
isotropic  and  linearly  elastic. 

The  earlier  fracture  algorithms  are  very  general  with  re- 
gard  to  describing  the  strength  distribution  of  a  subvolume  of 
material  under  arbitrary  temperature  and  state  of  stress.  In 
fact  they  are  so  general  as  to  render  them  virtually  impossible 
to  directly  apply  in  all  but  the  simplest  situations.  The 
governing  generalized  stress  has  to  be  determined  and  then  the 
strength  distribution  must  be  determined  for  each  temperature 
with  no  requirement  that  the  distributions  be  of  any  particular 
form  or  that  the  same  form  be  preserved  throughout  the  tempera¬ 
ture  range  of  interest. 


In  this  formulation  of  the  fracture  algorithm,  it  is  as¬ 
sumed  that  the  behavior  of  the  material  under  a  general  two- 
dimensional  state  of  stress  can  be  adequately  described  by  the 
approximate  biaxial  statistical  fracture  theory  developed  in 
Appendix  I.  In  this  theory,  the  probability  of  failure  of  a 
subvolurae  with  volume  AV  under  a  biaxial  state  of  stress 


o 


1' 


is  given  by 


exp 


AV 

V 


[f(aO 


f(cf2)j  ) 

i 


where 


f(a) 

f(a> 


0-C ( 


,  m 


u 


for  a  ■>  a, 


u 


for  o*au 


(1) 


(2) 


This  formulation  has  the  desirable  property  that,  for  the  case 
of  uniaxial  tension  (o^  *  0),  it  reduces  to  the  familiar  Weibuil 
form.  Consequently  all  the  procedures  developed  for  determining 
the  Weibuil  parameters  can  be  utilized  -  providing  that  all  the 
strength  distributions  are  of  the  Weibuil  form  with  perhaps  the 
parameters  varying  with  temperature.  It  should  be  noted  that 
the  assumption  expressed  by  Equation  (1)  does  not  automatically 
imply  that  the  material  is  a  pure  series  material. 


When  a  material  is  not  known  to  be  a  pure  series  material 
some  additional  care  must  be  exercised  regarding  the  size  of 
the  subvolumes  into  which  a  structure  is  divided.  The  reason 
for  this  is  that  if  a  material  is  not  pure  series,  then  a  non¬ 
conservative  step  is  performed  if  the  volume  of  the  subvolumes 
iV  are  made  smaller  than  the  tensile  strength  distribution 
specimen  volume  v  in  that  the  strength  of  the  subvolume  is 
overestimated.  In  order  to  avoid  this  problem  we  suggest  that 
all  subvolumes  be  made  equal  to  or  larger  than  the  tensile 
specimen  volume,  i.e.,  V>  v  for  all  subvolumes. 

3.  FRACTURE  ALGORITHM 

The  following  statement  of  the  fracture  algorithm,  along 
with  the  aforementioned  assumptions  has  ueen  devised  specific¬ 
ally  for  the  analysis  of  the  structures  tested  on  this  pro¬ 
gram. 

(1)  At  each  temperature  obtain  the  strength 
distribution  curve  F~(a)  for  a  tension 
specimen  of  volume  v.  Fit  the  distribution 
curve  with  the  form 

,  )  v 

L-exn  < 


F-(tf)  »  1-exp 


l 


v 

V 


0., 


for  oscr. 


u 


(3) 


m 


In  this  manner  the  parameters  are  determined 
as  functions  of  temperature:  m  w  m(T),  c?u  «  ou(T) 
and  o0  «  o0(T), 

For  each  mechanical  loading  and  thermal  loading 
(temperature  distribution)  determine  the  stress 
distribution  throughout  the  structure: 
a  -  s  a{  (x,y,2,t),  i  =>  1,2 

A  *. 
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(3) 


(4) 


Bi.(t)  = 


(5) 


(6) 


Divide  the  structure  into  n  convenient  sub¬ 
volumes  V, ,  V V  .  No  subvolume  should  be 
smaller  than  the  gage  volume  of  the  tension 
specimen;  v,  j  =  l,2,...n.  Subvolumes 
should  be  selected  with  approximately  homo¬ 
geneous  stress  states. 

For  each  value  of  time  t  determine  the 
"worst"  risk  of  rupture  for  each  subvolume  Vj  : 


V.  (  r 

v  f  m(T),  au( T),  aQ 


where  f 


I"  ai(x,y,z,t)  -  au(T)] 


m(T) 


for 


°fau 


,1=1,2 

(4  ) 

(5  ) 


=  0  for  a^<  ou 


and  where  T  a  T(x,y,z,t)  and  Dj  is  the  region 
in  space  occupied  by  . 

Determine  the  maximum  value  of  each  B^(t)  in 
the  interval  O-st^T  : 

bV  j  (T )  “  max  B^ .  (t)  (  6  ) 

*3  Q<  3 


The  probability  that  the  entire  structure  will 
survive  the  entire  environmental  history  up 
to  t  «  t  is  given  by 


where  the  term 

exp  |  -  [  By(t)  +  B2j  (t)]  \ 

^  ch 

represents  the  reliability  of  the  j 

subvolume, 1-Fj . 
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4.  EXAMPLE 

Consider  the  problem  of  calculating  the  reliability  of  the 
circular  plate  shown  in  Figure  1  .  Although  the  plate  is  as¬ 

sumed  to  be  at  uniform  temperature  and  the  mechanical  loading 
does  not  vary  with  time,  the  problem  should  still  illustrate 
several  of  the  facets  in  the  application  of  the  fracture 
algorithm.  Following  the  algorithm  we  proceed  as  follows. 

(1)  For  the  purposes  of  this  example,  the  tensile 
strength  distribution  is  assumed  to  be  of  the 
Weibull  type  with  parameters:  m=3.0,  c?u  =  4,130  psi 
and  aQ  =  5,610  psi  and  the  tension  gage  volume  is 
v  =  0.0125  in? 


(2)  Assuming  a  linearly  elastic,  homogeneous 
and  isotropic  material,  the  radial  and 
circumferential  stresses  are  found  from 
elasticity  theory  to  be 


o  °  S. 
r  i 


m-j  1 1 ri  i2 

”3 - r  •“/  +  solT - 2  i“ 

L“-o  •  ri  J  h0-rlJ 


2  2  ' 
r‘  +  r4  i/r 


ro  *  ri  J 


!1U2  + 

v  r  i 


2,2 
r  +  r . 
x 


-  r . 2  J  1  r  ' 


For  S. 


4,956  psi 


«  4,460  psi 

r  “  4.0  in. 
o 

r^  *=  1.0  in. 

h  -  1.0  in. 

these  expressions  reduce  to 

i  nr,  0.495 
O  =  4.99  -  - ij— 

• 


Oq  "  4.99  + 


0.495 
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Figure  1  Circular  Plate  with  Central  Hole  Loaded 

by  Both  Inside  and  Outside  Tensile  Stresses 


(3)  The  plate  is  divided  into  five  concentric 
rings  with  the  dimensions  indicated  in 
Table  I  . 

(a)  Note  that  each  ring  volume  is  greater 

_  3 

than  the  gage  volume  v  =  0.125  in. 

(b)  The  stress  state  in  each  ring  becomes 
more  homogeneous  when  the  number  of 
rings  is  increased. 

(4)  Due  to  the  fact  that  the  temperature  distribu¬ 
tion  is  uniform,  the  "worst"  risk  of  rupture 
is  computed  using  the  maximum  stresses  in  the 
subvolume.  The  largest  radial  stress  in  the 

j  L  ring  occurs  at  its  outside  radius  bj 


arj  =4.99  - 


0.495 


The  largest  circumferential  stress  is  found 
at  the  inside  radius  a. 


"  4'99  + 


0.495 


Thus  the  components  of  the  "risk  of  rupture" 
th 

for  the  j  ring  are  given  by 


*  •»' 
J8i  ’  qu 

% 


The  values  of  the  stresses  and  "risks  of 
rupture"  are  tabulated  in  Table  I  . 
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(5)  Due  to  the  fact  that  there  is  no  variation 
with  time,  we  have 


(6) 

l-Fj  -  exp  [-  b”j  -  B~j  05) 

and  tabulating  the  results  in  Table  I  ,  the 
reliability  of  the  plate  is  computed  to  be 
5 

1-F  *  H  (1-F. ) 
j-1  J 

»  (0.94148)  (0.94996)  (0.94090) 

(0.93674)  (0.85015)  (16) 

-  0.670 

If  the  material  in  ! he  plate  is  known  to  be  a  series 
material,  we  can  drop  the  restriction  that  v  .  This 

enables  us  to  use  infinitesimal  rings  which  leads  to  a  reliability 
prediction  of  0.696.  Thus,  in  this  example,  the  partitioning 
of  the  plate  into  only  five  unit  volumes  results  in  a  fair 
estimate  of  the  total  reliability. 


ic 

B  .  =  B  . 

rj  rj 


(14) 


Br 


=  B, 


9j 

Computing  the  ring  reliabilities 
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SECTION  III 


EXPERIMENTAL  PROGRAM 

1.  BEAMS:  THERMAL  AND  MECHANICAL  LOADING 
a.  General  Discussion 


The  specimen  selected  for  f L  xural  testing  under  thermal 
and  mechanical  loading  was  a  beam  1/2  in.  wide,  1  in.  deep, 
and  10  in.  long  -  made  of  Wesgo  A1  995  material.  The  test  mode 
was  four  point  bending  with  a  gage  span  of  4  in.  Because  of 
the  additional  condition  of  thermal  loading  it  was  necessary 
to  construct  a  loading  fixture  which  would  bridge  the  heating 
apoaratus  (Figure  2).  Further,  loading  was  done  through  a 
point  contact  and  roller  scheme  to  compensate  for  any  initial 
misalignment  in  the  specimen  (see  detail  in  upper  corners  of 
Figure  2) . 

The  thermal  loading  was  obtained  through  the  use  of  two 
strip  heaters  (R  I  Controls  model  number  5305- 5A) .  The  units 
consisted  of  parabolic  reflectors  which  were  focused  to  concen¬ 
trate  the  heat  flux  onto  the  tension  side  of  the  beam.  This 
arrangement  necessitated  protecting  the  quarts  heating  element 
from  broken  beam  debris  with  a  piece  of  heavy  screen.  Antici¬ 
pating  some  change  in  heat  flux,  all  temperature  distributions 
were  obtained  with  this  screen  in  place.  When  preparations  were 
complete,  several  beams  were  broken  without  the  thermal  loading 
as  part  of  a  general  check  on  the  system,  and  the  flexural 
strengths  were  found  to  be  in  the  expected  range, 

b.  Temperature  Distribution 

The  temperature  distribution  as  a  function  of  time  was  ob¬ 
tained  with  the  use  of  rapid  response  iron  constantan  thermo¬ 
couples  and  an  electronic  multiple  channel  millivolt  recording 
instrument,  (each  channel  prints  every  12  seconds).  Six  1/16-in. 
diameter  holes  were  drilled  into  the  depth  of  the  beam  from  the 
top  side  at  varying  depths  of  0.179,  0.314,  0.451,  0.552,  0.727, 
and  0.363  in.  (Figure  3).  The  temperatures  were  monitored  with 
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the  beam  in  place  in  the  test  fixture  duplicating  the  conditions 
which  would  be  encountered  during  the  actual  testing  of  the  speci¬ 
mens.  A  refractory  brick  was  notched  to  the  width  of  the  beam  to 
provide  the  insulation  and  shielding  required  to  discourage  all  but 
a  pure  one-dimensional  heat  flow  through  the  beam  depth  (Figure  4) . 
The  recording  instrument  was  calibrated  and  the  thermocouples 
were  checked  in  a  water  bath  at  boiling  point. 

Initial  trials  revealed  that  the  A^Q-j  material  exhibited 
some  transparency  to  thermal  radiation  which  distorted  the 
reading  of  those  thermocouples  nearest  the  hot  interface.  This 
problem  was  dealt  with  by  coating  this  surface  with  a  thin  even 
layer  of  nickel  silicate.  Four  sets  of  data,  each  determining 
a  temperature  distribution,  were  obtained.  Each  set  was  checked 
against  the  other  to  ascertain  reproducibility.  The  temperature 
distribution  curves  are  shown  in  Figure  5. 

c.  Test  Procedure 

The  test  procedure  consisted  of  applying  an  initial  load  of 
50  lb  and  energizing  the  strip  heaters  at  a  controlled  level  of 
190  volts  for  8  min.  After  8  min.  the  mechanical  load  was  in¬ 
creased  at  an  approximate  rate  of  500  lb  per  min.  until  fracture 
occurred.  The  time  to  failure  is  taken  to  be  8  min.  plus  the 
duration  of  increasing  load.  The  initial  50  lb  loading  served 
to  preserve  initial  alignment  and  eliminate  backlash  during 
change  from  thermal  to  thermal-mechanical  loading.  A  total  of 
48  beams  were  tested  in  this  manner.  The  test  machine  was  equipped 
with  a  load  pacer  device  which  enabled  the  operator  to  approximate 
a  predetermined  load  rate.  To  compute  the  actual  load  rate,  the 
time  from  the  start  of  the  load  increase  to  final  fracture  was 
recorded  for  each  test  and  loading  rates  were  determined  by  divid¬ 
ing  the  final  load  minus  50  lb  by  the  duration  of  the  increase. 
Consequently,  there  was  a  loading  rate  associated  with  each  test 
which  deviated  somewhat  from  the  target  number  of  500  lb  per  min. 


From  this  data  an  estimate  of  the  average  loading  rate  for  all 
tests  was  determined  as  492.9  lb  per  min.  This  number  was  in  turn 
used  to  determine  a  corrected  time  to  failure. 

The  statistical  significance  of  these  operations  is  based  on 
our  hypothesis  that  the  mechanical  and  thermal  load  histories  are 
the  same  for  all  specimens;  this  requires  that  the  mechanical 
loading  rate  be  the  same  for  each  beam.  The  actual  load  rates 
are  not  constant  and  do  not,  therefore,  exactly  conform  to  our 
hypothesis.  To  account  for  the  small  variations  in  loading  rate, 
we  used  the  average  load  rate  to  compute  a  corrected  time  to  fail¬ 
ure  which  would  reflect  the  actual  failure  load.  This  correction 
is  possible  only  because  the  thermal  stresses  in  the  beam  do  not 
significantly  contribute  to  its  probability  of  failure.  (See 
Section  IV-4) . 

The  beam  fracture  data  is  presented  in  Table  II.  Column  1 
tabulates  in  ordered  form  the  fracture  time  recorded  for  each 
test.  The  loading  rate  for  each  test  is  listed  in  column  5  and 
the  average  of  column  5  which  is  the  average  loading  rate  is  de¬ 
termined  to  be  492.9  lb  per  min.  The  corrected  time  to  failure 
is  found  by  dividing  column  4  by  492.9  and  adding  S.O  min. 

Table  III  presents  the  corrected  and  ordered  data  and  Figure  6 
illustrate  the  resulting  cumulative  distribution. 

2.  PLATES:  THERMAL  LOADING 

a.  General  Discussions 

The  plate  specimen  waj  made  of  Wesgo  A1  995.  The  configura¬ 
tion  was  circular  with  a  concentric  1-in.  diameter  hole.  The 
overall  diameter  was  6  in.,  and  the  thickness  dimension  was  1/4  in. 

The  thermal  loading  produced  by  means  of  a  3/4- in.  diam¬ 
eter  .silicon  carbide  heating  element  (Globar)  positioned  perpen¬ 
dicular  to  the  plate  through  its  center  hole.  A  purely  two- 
dimensional  heat  flow  was  encouraged  by  insulating  both  the  top 
and  bottom  of  the  plate  so  that  the  principal  heat  loss  would  be 
through  the  outer  place  edge.  Pictures  of  this  setup  are  shown 
in  Figures  7  and  8. 
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Table  II 


TEST  DATA  -  A£203  BEAM 


1 

2 

3 

4 

5 

6 

TIME  TO 

FAILURE 

TIME  TO 

failure 

CORRECTED 

TIHE  TO 

FAILURE 

LOAD 

FAILURE 

LOAD 

LOADING 

FAILURE 

(MIN.) 

(LBS.) 

MINUS 

8  MIN. 

MINUS 

50  LB. 

RATE 

COL  4  /  COL  3 

MINUS  8  MIN. 
COL  4/  492.9 

..  JOB, 

1680 

3.283 

1630 

496 

3.307 

. 11.283.  .. 

.  1 670  ___ 

.  _  3,283. _ 

_ 1620 _ 

493 

3.29 

.  11.384. 

1755 

. 3. 384 

1705 

_ 504 _ 

3.459 

_ U.8.,  .  . 

...  1.820 _ 

3.6 

1770 

491 

3.692 

_  l!,7$  . 

1940 _ 

3.75 

1890 

504 

3.834 

_ UJ!8S.. 

_ 1320 _ 

3.784 

1070 

494 

3.780 

11.95 

2000 

1950 

3.956 

11.95 

2010  ~ 

i960— _ 

3.976 

11.95 

2040 

_ 3.95_ 
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Table  III 


AI0O3  BEAM  DATA  -  CORRECTED  TIME  TO  FAILURE 
AND  ASSOCIATED  PROBABILITY 
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Figure  6  Cumulative  Distribution  of  BEAM  DATA 
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Figure  8  Cross  Section  Showing  Detail  of  Various 
Test  Components 


The  radial  loading  fixture  which  appears  in  all  of  these  photo- 
graphs  is  present  only  for  the  mechanical  aspect  of  the  loading 
which  is  disucssed  under  part  3  of  this  section.  The  power  input 
to  the  heating  element  was  measured  and  controlled  by  the  use  of 
a  watt  meter  and  rheostat.  The  test  was  completed  when  catas¬ 
trophic  failure  of  the  disk  occurred.  The  time  to  failure  in 
each  test  was  the  interval  between  the  instant  of  power  applica¬ 
tion  and  the  instant  of  failure.  The  gross  heat  flux  incident 
upon  the  disk  as  computed  using  the  temperature  distribution  and 
material  properties  will  be  presented  later  in  Section  IV-3. 

b.  Temperature  Distribution 

As  in  the  case  of  the  beam  it  was  necessary  to  obtain  a  temp¬ 
erature  distribution  throughout  the  disc  as  a  function  of  time 
and  location.  Because  of  the  axisymmetric  character  of  the  ther¬ 
mal  loading,  the  location  is  specified  by  radius  only.  A  very 
satisfactory  way  of  monitoring  temperatures  through  the  disk  con¬ 
sisted  of  implanting  thermocouples  at  strategic  locations  on  the 
plate  (Figures  9  and  10) .  This  precluded  the  possibility  of  in¬ 
strumenting  each  plate  tested.  Hence  it  was  necessary  to  obtain 
a  definitive  set  of  measurements  from  one  fully  instrumented  disk 
to  infer  the  temperature  distribution  which  existed  during  all 
successive  tests.  As  in  the  case  of  the  beam,  some  transparency 
near  the  inside  hole  was  in  evidence  and  a  graphite  coating  was 
applied  to  the  inside  edge  of  ail  plates  to  minimize  "see  through." 

It  should  be  recalled  that  failure  is  catastrophic  (Figure  11) 
in  these  tests  and  that  the  time  to  failure  is  variable.  Hence, 
the  problem  arises  of  obtaining  a  representative  set  of  temperature 
data  for  a  duration  of  time  greater  than  the  greatest  time  to  fail¬ 
ure  expected  for  ail  tests  without  having  the  instrumented  disk 
fail  during  this  period.  The  only  solution  that  presented  itself 
was  to  "prefail"  the  instrumented  disk.  Two  pieces  were  formed 
by  cracks  along  the  radial  direction.  See  Figures  12  and  13  for 
explanatory  diagrams  showing  location  of  the  cracks  and  thermo¬ 
couple  implants.  Several  preliminary  tests  on  plate  specimens 
gave  an  indication  of  what  to  expect  as  a  time  to  failure  duration. 


27 


Figure  9  General  View  of  the  Instrumented  Plate  in  Test  Mode 


Detailed  View  of  Instrumented  Plate  with  Insulation 
Removed  to  Show  Thermocouple  Implants 


Figure  11  Typical  Failure  Mode  After  Application  of  Thermal  Load 


of  Thermocouples 


Using  this  as  a  guide,  the  instrumented  disk  was  monitored  for  a 
time  period  almost  twice  that  expected.  Finally,  many  tests  were 
run  with  this  instrumented  disk  under  a  variety  of  circumstances 
which  served  two  primary  interests.  First,  that  of  obtaining  a 
temperature  distribution  which,  with  a  degree  of  confidence,  will 
accurately  represent  the  actual  distribution  in  each  test  speci¬ 
men.  Secondly,  planned  differences  such  as  dismantling  setup  and 
reassembling,  immediate  reruns  without  any  disturbance  and  rota¬ 
tional  reorientation  of  the  disk  with  respect  to  the  heating  element 
gave  the  technicians  experience  which  was  a  valuable  asset  in  con¬ 
ducting  the  unmonitored  disk  experiments. 

The  temperature  data  was  obtained  as  a  millivolt  reading  at 
a  particular  time  and  position.  This  information  has  been  con¬ 
verted  to  temperature- time-position  measurements  and  replot. rad 
in  a  manner  more  useful  in  the  program.  The  results  of  these 
temperature  distribution  tests  are  shown  in  Figure  14. 

c.  Fracture  Testing 

Results  of  these  tests  are  shown  in  Table  IV  and  the  result¬ 
ing  cumulative  distribution  is  illustrated  in  Figure  15. 

3.  PLATES;  THERMAL  AND  MECHANICAL  LOADING 

a.  General  Discussion 

When  we  formulated  our  original  program  strategy,  it  was  con¬ 
templated  that  the  disk  element  would  be  subjected  to  an  axisym- 
metric  thermal  loading  and  a  mechanical  radial  tensile  loading. 
Unfortunately,  the  required  magnitude  for  the  tensile  loading 
could  not  be  achieved  with  the  gripping  scheme  that  we  visualized. 
For  this  reason,  it  seemed  expedient  to  consider  a  radial  com¬ 
pressive  loading  which  would  not  present  special  gripping  problems. 

Using  the  radial  load  fixture  shown  in  Figures  16  and  17, 
compressive  loads  were  applied  to  the  disk  through  various  types 
of  bearing  devices.  We  attempted  to  approach,  as  close  as  pos¬ 
sible,  a  uniform  peripheral  loading.  Our  efforts  were  simply  not 
good  enough  as  shown  by  the  photoelastic  results  in  Figure  18. 
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Table  IV 


AijO,  DISK  DATA  -  TIME  TO  FAILURE  AND  ASSOCIATED  PROBABILITY 

(Thermal  Loading  Only) 
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Figure  15  Cumulative  Distribution  of  A_0~O~  Disk  Data  (Thermal  Load  Only) 


Plate  During  Application  of  Both  Mechanical  and 
Thermal  Load 


Figure  17  Mechanically  Loaded  Plate  (Coinpres 


Here,  we  observe  not  only  that  the  stress  distribution  is  not  axi- 
symmetric;  but,  that  significant  tensile  stresses  are  present.  As 
we  shall  see,  these  tensile  stresses  cause  the  disk  to  fail  sooner 
than  the  disks  without  this  compressive  loading.  This  illustrates 
quite  clearly  that  a  compressive  prestress  will  not  necessarily 
strengthen  a  brittle  element.  We  hasten  to  point  out  that  a  uni¬ 
form  compressive  prestress  would  not  have  introduced  tensile  stres¬ 
ses  and  would  have  resulted  in  an  increase  in  the  plate's  integrity. 

b.  Temperature  Distribution 

The  temperature  distribution  for  the  disk  with  compressive 
grips  in  place,  was  obtained  in  the  same  manner  previously  de¬ 
scribed  in  part  2b  of  this  section.  The  temperature  distribution 
curves  are  shown  in  Figure  19.  Because  90  percent  of  the  periphery 
was  used  to  approach  a  uniform  mechanical  loading,  the  effective 
heat  sink  was  appreciable.  Comparison  with  the  previously  obtain¬ 
ed  temperature  distribution  shows  a  strong  shift  to  lower  tempera¬ 
tures  for  the  same  time  intervals,  as  expected. 

c.  Fracture  Testing 

The  disks  were  placed  in  the  loading  fixture  as  shown  in  Fig¬ 
ure  16  and  subjected  to  a  mechanical  loading  and  thermal  loading, 
the  latter  being  superposed  on  the  mechanical  loading  when  it 
reached  a  prescribed  level.  The  time  to  failure  was  recorded  for 
each  test  and  is  given  in  Table  No.  V.  Failure  was  defined  by  the 
first  audible  sound  of  cracking.  The  cumulative  probability  dis¬ 
tribution  is  shown  in  Figure  20. 


’able  V 


3K  DATA  -  TIME  TO  FAILURE  AND  ASSOCIATED 


OBABILITY  (Mechanical  and  Thermal  Loading) 
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Cumulative  Distribution  of  A£00„  DISK  DATA  (Thermal  and  Mechanical  Loading) 


SECTION  IV 
RESPONSE  PREDICTIONS 


1.  APPROACH 

The  central  question  to  which  this  research  effort  addresses 
itself  is:  "Can  the  statistical  analysis  algorithm  be  successfully 
applied  in  the  case  of  brittle  structures  exposed  to  thermal  and 
mechanical  loadings  and,  if  not,  why  not?"  Specifically,  on  this 
program  we  seek  to  explain  the  results  of  the  beam  and  disk  ex¬ 
periments  described  in  Section  III  with  the  aid  of  the  fracture 
algorithm  presented  in  Section  II. 

There  are  many  possible  avenues  of  approach  on  this  question. 
The  straightforward  approach  is  the  most  logical  and  involves 
using  experimentally  obtained  basic  material  property  data  in 
conjunction  with  the  fracture  algorithm  to  attempt  to  "predict" 
the  results  of  the  beam  and  disk  experiments.  The  difficulty 
with  this  approach  is  that  the  tensile  strength  distribution  for 
alumina,  at  room  temperature  and  especially  at  elevated  tempera¬ 
tures,  has  not  yet  been  satisfactorily  defined.  Each  available 
strength  distribution  is  probably  an  underestimate  of  the  true 
distribution  for  alumina.  Even  if  alumina  were  a  pure  series 
material  of  the  Weibull  type  and  the  algorithm  was  performed 
using  infinitesimal  subvolumes,  the  use  of  the  available  strength 
distributions  should  result  in  conservative  and  not  exact  results. 
Thus,  the  achievement  of  conservative  predictions  is  not  a  very 
sensitive  test  of  the  applicability  of  fracture  algorithm  - 
especially  if  the  material  is  close  to  being  of  the  pure  series 
type  as  in  fact  we  hoped  alumina  would  be. 

The  approach  which  was  selected  begins  by  making  an  additional 
assumption  that  the  tensile  strength  distribution  for  alumina  is 
independent  of  temperature  in  the  range  from  R.T,  to  15Q0°F.  At 
moderate  temperatures  some  investigators  (References  3, 4, 5, 6  and  7) 
do  measure  a  slight  degradation  of  strength  with  increasing 
temperature.  However,  our  experience  with  elevated  temperature 
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testing  suggests  that  the  experiments  cannot  be  performed  with 
the  same  precision  as  they  can  at  room  temperature.  Thus,  since 
a  measured  degradation  of  strength  at  moderate  temperatures  may 
be  explained  by  experimental  limitations,  the  assumption  of  the 
strength  distribution  being  temperature  independent  may  not  be 
unreasonable. 

Since  the  strength  distribution  is  assumed  to  be  of  the 
Weibull  type,  the  assumption  of  temperature  independence  requires 
that  the  parameters  m,  a  ,  oq  are  independent  of  temperature. 

In  this  approach  these  parameters  are  left  as  open  parameters. 

A  trial  and  error  procedure  employing  the  fracture  algorithm  in 
conjunction  with  the  stress  analysis  is  used  to  attempt  to  simul¬ 
taneously  "fit"  both  the  beam  and  disk  experiments  using  a  single 
set  of  values  for  the  parameters.  Using  this  method,  a  deficiency 
in  the  algorithm,  experimental  procedure,  or  the  temperature  in¬ 
dependent  strength  assumption  may  be  observed  in  two  ways.  First, 
it  is  entirely  possible  that  the  beam  and  disk  experiments  cannot 
be  "fitted"  satisfactorily  using  a  single  set  of  parameters.  Sec¬ 
ond,  assuming  that  it  is  possible  to  fit  both  the  beam  and  disk 
experiments,  it  is  very  likely  that  the  resulting  strength  distri¬ 
bution  may  be  unreasonably  different  from  the  data  obtained  from 
the  best  strength  tests . 

2.  MATERIAL  PROPERTIES 

The  approach  just  described  was  successfully  employed.  Using 

the  single  set  of  Weibull  parameters  m  a  3.3,  a  a  19,700  psi, 

“  3 

crQ  «  7,000  psi  with  a  tension  gage  volume  of  v  «  0.0982  in., 

both  the  beam  and  disk  experiments  were  satisfactorily  "fitted." 
Using  these  values  the  strength  distribution  is  plotted  in  Fig¬ 
ure  21.  The  mean  strength  can  be  computed  to  be  approximately 
32,400  psi.  This  value  tends  to  be  somewhat  higher  than  most  avail¬ 
able  tension  data  and  begins  to  approach  the  flexural  strength  data. 
However,  very  recent  gas-bearing  data  by  Pears  and  Starrett  (Ref¬ 
erence  8)  achieves  some  high  mean  strengths  in  tension  for  alumina. 
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When  the  comparison  is  made  between  this  excellent  data  for  a  gage 
volume  of  0.031  in.  to  the  derived  Weibull  type  strength  dis¬ 
tribution  for  the  same  volume,  Figure  22  ,  the  agreement  is 
seen  to  be  quite  remarkable. 

For  the  remainder  of  this  report  the  tensile  strength  dis¬ 
tribution,  in  the  range  from  R.T.  to  1500°  F,  shall  be  assumed 
to  be  given  by  the  Weibull  formulation  with  parameters  m  =  3.3, 
ou  =  19,700  psi  and  aQ  =  7,000  psi.  Furthermore,  for  the  pur¬ 

pose  of  illustrating  the  general  application  of  the  fracture 
algorithm,  the  distribution  curve  shall  be  pretended  to  come 
from  tests  on  a  tensile  specimen  of  gage  volume  v  =  0.0982  in. 

The  other  material  properties  entering  into  the  stress 
analysis  and  thermal  analysis  have  been  obtained  from  the  litera¬ 
ture.  Thermal  strain  as  a  function  of  temperature  after  Goldsmith, 
Hirschhora  and  Waterman  (Reference  9)  is  shown  in  Figure  23.  Mod¬ 
ulus  of  elasticity  as  a  function  of  temperature  after  Dally  (Ref¬ 
erence  10)  is  illustrated  in  Figure  24.  Poisson's  ratio  as  a 
function  of  temperature  after  Coble  and  Kingery  (Reference  11)  is 
shown  in  Figure  25.  Finally  the  thermal  conductivity  as  a  function 
of  temperature  after  Goldsmith,  Hirschhom  and  Waterman  (Reference  9) 
is  illustrated  in  Figure  26. 

3.  THERMAL  ANALYSIS 

Assuming  black  body  radiation,  the  heat  flux  Q  impingent 
upon  the  surface  of  the  beam  or  disk  may  be  expressed  as 

Q  “  -  k  +  o  T4  (17) 

^T 

where  k  is  the  thermal  conductivity  of  alumina,  — is  the 
normal  derivative  of  the  temperature  distribution  into  the 
body  evaluated  on  the  surface,  a  is  the  Stefan-Boltzman  con¬ 
stant  and  T  is  the  absolute  temperature  on  the  surface.  Using 
the  temperature  distribution  curves  for  the  beam  and  disk  from 
Section  II,  Figures  5  and  13  ,  and  Figure  26  for  the  thermal 
conductivity,  the  heat  flux  Q  can  be  computed  using  Equation  (17). 
The  results  of  these  computations  are  plotted  in  Figure  27. 


Figure  22  Comparison  of  Tensile  Strength  Cumulative  Distribution  for  Alumina  to 
Gas-Bearing  Data  of  C.  D.  Pears 


Thermal  Strain  of  Alumina  as  a  Function  of  Temperature  (After  Goldsmith 
Hirschhom  and  Waterman) 
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Figure  27  Gross  Heat  Flux  Incident  Upon  Beam  and  Disk 
as  a  Function  of  Time 
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4.  BEAM  STRESS  ANALYSIS 


.  The  beam  prototype  structure  is  loaded,  as  shown  in  Figure  28, 
by  terminal  couples  M(t)  and  by  the  temperature  distribution  T(z,t) 
induced  by  the  heat  flux  Q(t)  impinging  on  the  bottom  surface.  The 
general  solution  of  the  thermal  stresses  in  a  beam  with  temperature 
dependent  mechanical  and  thermal  properties  using  strength-of-ma- 
terials  type  assumptions  is  well  known  (Reference  12) .  By  replacing 
the  term  aT  everywhere  by  the  thermal  strain,  and  introducing 
the  mechanically  induced  bending  moment  distribution  M(t) ,  the  stress 
distribution  a.  (z,t)  is  found  to  be  given  by 


pT(t)  [ttj  (t)  +  M(0  > 


a 


e^Ly  d  A 


It  is  observed  that,  in  this  formulation  with  e^,  replacing 
aT,  the  dependence  of  the  solution  upon  the  temperature  dis¬ 
tribution  is  now  entirely  implicit. 

Due  to  the  fact  that  the  temperature  distribution  is  in¬ 
dependent  of  y  the  following  simplifications  result: 


(22) 


Mj  (t)  =  0 

2  d 


PT(t)  =  b 


SipE  d  z 


Mj  (t)  =  b  eTE  z  d  z 

y  J 

-d 

3 

Noting  th-t  A  =  2bd  and  I  =  (2/3)bd  ,  the  expression  for  Che 
stress  distribution  becomes 

rd  A 

3z 

lT“  '2d 


a  (z,t)  =  -  e^E  +  1 

x~  2bdJ 


£p,Ed2  + 


2d 


£tE  z  d  z  (23) 


-d 


Thus,  given  £p  ®  £j(T),  E  =»  E(T),  T(z,t)  and  M(t) ,  Equation  (23) 
will  yield  the  stresses  in  the  beam. 

The  solution  of  Equation  (23)  has  been  accomplished  using 
a  digital  computer  (Appendix  II).  Using  the  temperature  distribu¬ 
tion  curves,  Figure  5  ,  the  material  property  curves,  Figures  23 
and  24,  and  the  relationship 


M(t)  °»  59.4  in-lbs  for  Oat  *8.0  min 

«  59.4  +  535.3  (t-3.0)  in.  lbs  for  t  8.0  min 


(24) 


Che  stress  distribution  is  found  to  be  as  shown  in  Figure  29. 

It  is  observed  that  the  stresses  are  small  until  after  the  me¬ 
chanical  load  begins  to  increase  at  the  8.0  min  point.  Also, 
another  consequence  of  the  dominating  mechanical  loading  is  that 
the  stress  distribution  is  v ~ry  nearly  linear  after  the  9.0  min 
point. 

Figure  30  illustrates  the  maximum  tensile  stress  distribu¬ 
tion.  It  is  interesting  to  note  Char,  using  the  mean  tensile 
strength  of  approximately  32,-00  psi,  the  deterministic  maximum 
stress  theory  predicts  failure  at  about  12.6  min.  From  the  ex¬ 
perimental  data,  Figure  6,  the  earliest  failure  occurs  at  about 
11.3  min  and  the  median  time  to  failure  is  about  12.3  min. 
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5 .  DISK  STRESS  ANALYSIS 


The  disk  prototype  structure  is  loaded,  as  shown  in  Figure  31, 
by  a  tensile  radial  stress  S  on  the  outside  surface  and  by  the 
temperature  distribution  T(r,t)  produced  by  the  heat  flux  Q(t)  on 
the  inside  surface.  The  plane  stress  solution  in  the  case  where 
the  mechanical  and  thermal  properties  of  the  material  are  independ¬ 
ent  of  temperature  is  well  known  (Reference  12) .  One  method  of 
accounting  for  temperature  dependent  properties  is  to  represent 
the  disk  as  a  composite  structure  of  concentric  rings  where  each 
ring  is  thin  enough  so  that  the  material  properties  within  the 
ring  can  be  assumed  to  be  constant. 

The  analysis  of  the  concentric  rings  can  be  readily  accom¬ 
plished  using  the  Method  of  Initial  Parameters  or  Transfer  Matrix 
Method  (Reference  13)  due  to  the  fact  that  the  geometric  and  equi¬ 
librium  conditions  which  must  be  satisfied  at  each  interface  be¬ 
tween  rings  are  automatically  fulfilled.  In  addition  to  those 
parameters  defined  in  Figure  31  ,  the  parameters  and  variables 
entering  into  the  problem  are: 

Em  «  thermal  strain 


T(r.t) 

u(r,t) 


ar(r,t) 


oQ(r,t) 


modulus  of  elasticity 

Poisson’s  ratio 

temperature  distribution 

radial  displacement  distribution 

radial  stress  distribution 

circumferential  stress  distribution 


Using  a  matrix  form  of  this  method,  the  stress  distribution  in 
the  disk  is  given  by 


a  (r,t)  ® 


Nr(r,t) 


oQ(r,t)  » 


T - 

N0(r,t) 
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Figure  31  Concentric  Ring  Geometry  for  Disk  Analysis 


where 
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In  general  Equation  (26  )  leads  to  two  nontrivial  relation¬ 
ships  which  may  be  expressed  as 


u(r,t)  ■  uQ(t)  Fj^Cr.t)  +  NQ(t)  F2(r,t)  +  F3(r,t) 

(2  2) 

Nr(r,t)  -  uQ(t)  GL(r,t)  +  NQ(t)  G2(r,t)  +€3(r,t) 


where  UQ(t)  =  u(rQ,t)  and  NQ(t)  =  Nr(rQ,t)  and  where  F^(r,t) 
and  G^(r,t),  i  =  1,2,3,  represent  the  cumbersome  expressions 
that  could  be  obtained  from  Equation  (26) .  Now  since  the 
boundary  conditions  are 


Nr<r0>t)  =  0 
Nt(ri»t)  =  Sh 


(33) 


the  two  constants  of  integration  uQ(t)  and  NQ(t)  are  found  from 
Equation  (32)  to  be 


NQ(t) 


Sh  -  G3(r^,t) 

G^Cr^t) 


(34) 


The  system  of  equations,  Equations  (25)  through  (34),  has 
been  programmed  on  a  digital  computer  (Appendix  III) .  Using 
the  temperature  distribution  in  the  disk,  Figure  13,  and  the  curves 
for  the  material  properties,  Figures  23,  24  and  25,  the  stress 
distribution  has  been  computed  as  depicted  in  Figures  32  and  33. 

It  is  observed  that  the  radial  stress  is  always  compressive  and 
hence,  according  to  the  fracture  algorithm  does  not  contribute  to 
the  probability  of  failure.  The  circumferential  stress  distribu¬ 
tion  is  observed  to  be  compressive  near  the  inside  radius  and  ten¬ 
sile  over  most  of  the  remainder  of  the  disk. 

The  maximum  tensile  circumferential  stress  distribution  is 
illustrated  in  Figure  34.  Using  the  mean  tensile  strength  of 
approximately  32,400  psi,  the  deterministic  maximum  stress  theory 
predicts  that  no  disk  will  fail  up  to  the  end  of  the  test  at 
14.0  min.  From  the  experimental  data,  Figure  14,  the  earliest 
failure  occurs  at  just  under  7.0  min,  the  median  is  at  about 
9.6  min  and  all  specimens  have  failed  by  the  12.4  min  point. 
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6.  APPLICATION  OF  THE  FRACTURE  ALGORITHM 

Separate,  but  similar,  computer  programs  were  written  for 
the  analysis  of  the  beam  and  disk  prototype  structures  (see  Ap¬ 
pendices  II  and  III) .  Each  program  did  the  complete  analysis  of 
its  prototype  structure,  including  thermal  analysis,  stress 
analysis  and  the  statistical  analysis  employing  the  fracture 
algorithm. 

The  fracture  algorithm  described  in  Section  II  was  designed 
especially  for  the  analysis  of  the  beam  and  disk  experiments  and 
thus  no  major  modifications  were  necessary  in  incorporating  it 
into  the  computer  programs.  Steps  (3),  (4)  and  (5)  are  the  only 
steps  in  the  fracture  algorithm  that  have  not  already  been  dis¬ 
cussed  or  that  are  not  executable  in  their  given  form.  For  con¬ 
venience,  these  steps  shall  be  repeated  here  along  with  a  dis¬ 
cussion  of  how  the  steps  were  accomplished  in  the  analysis  of  the 
beam  and  the  disk. 

(3)  Divide  the  structure  into  n  convenient  subvolumes 
Vx,  ^2 , . . . .  No  subvolume  should  be  smaller  than 
the  gage  volume  of  the  tension  specimen: 

Vji  v,  j  »  l,2,...n.  Subvolumes  should  be  selected 
with  approximately  homogeneous  stress  states. 

The  volume  of  the  tension  specimen  is  0.0982  in.  The  depth 
of  the  beam  was  divided  into  20  equal  segments  thereby  creating 
identical  subvolumes  of  volume  0.100  in.  The  disk  was  subdivided 
into  69  ring  type  subvolumes,  each  with  a  volume  of  0.0996  in'*. 

By  selecting  the  smallest  volumes  possible  in  each  case,  the 
homogeneity  of  the  stress  state  in  each  subvolume  was  optimized. 

(4)  For  each  value  of  time  t  determine  the  "worst" 
risk  of  rupture  for  each  subvolume  V, : 


B. 


jU)  “  /max  f[<?t(x,y,z,t);  ra(T) , cru(T) , aQ(T)]  \  ,ial,2 
lDj  '  (4) 
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where 


at(x3y,z,t)  -  au(T) 
=  _  ^ 

=  0  for  at  2  au 


m(T) 


for  =.  a 


u 


(5) 


and  where  T  =  T(x,  y,z,t)  and 
space  occupied  by  . 


is  the  region  in 


In  the  computer  programs,  it  was  assumed  that  the  maximum 
value  of  f  always  occurred  on  the  boundary  of  the  subvolume  in 
both  the  case  of  the  beam  and  disk.  Regarding  the  beam  subvolume, 
the  value  of  f  was  computed  at  the  top  and  bottom  surfaces  and 
the  largest  was  selected  to  be  the  maximum  value  in  the  subvolume. 
In  the  case  of  the  disk,  f  was  computed  on  the  inside  radius  and 
outside  radius  with  largest  being  taken  for  the  maximum  in  the 
ring. 


Strictly  speaking,  the  maximum  value  of  f  will  not  always 
lie  on  the  boundary  of  the  subvolurae.  Now  Equation  (4)  is  itself 
a  conservative  statement  and  thus, when  the  assumption  that  f  Is 
maximum  on  the  boundary  is  valid,  the  results  are  still  conservative. 
However,  when  the  assumption  is  not  valid  the  results  are  not  con¬ 
servative.  Therefore  it  is  assumed  that,  by  making  the  subvolutnes 
small  enough,  more  often  than  not  the  maximum  values  will  be  at¬ 
tained  on  the  boundary  and  in  this  manner  the  cumulative  results 
for  the  risk  of  rupture  for  the  entire  structure  will  remain  con¬ 
servative  . 


(5) 


Determine  the  maximum  value  of  each  B^(t)  1°  Che 
interval  Of  t  ^  t 


* 


B,  »  (  c)  “  max 
O-t-T 


B^Ct) 


(6) 


In  the  computer  programs,  the  value  of  time  is  iterated  by 
finite  amounts  so  the  problem  of  determining  the  maximum  value  of 
•  B^j(t)  with  respect  to  time  is  similar  to  finding  the  maximum  of 
f  in  Dj  .  The  value  of  B^  is  computed  at  each  value  of  time 
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and  the  largest  value  is  selected  as  the  maximum  value.  Now  the 
true  maximum  value  may  have  occurred  at  some  intermediate  value 
of  time.  So ,  again,  an  appeal  is  made  to  small  increments  in 
order  to  minimize  the  error.  In  the  case  of  the  beam  and  disk, 
the  ' s  increased  fairly  monotonically  with  time  and  thus  the 
problem  was  alleviated  since  the  maximums  tended  to  occur  at  the 
current  value  of  time. 

The  remainder  of  tv  task  of  placing  the  analyses  on  the  com¬ 
puter  was  very  straightforward.  The  stress  analyses  were  pro¬ 
grammed  exactly  as  presented  earlier  in  this  section.  The  material 
properties  that  entered  into  the  analyses  were  incorporated  into 
the  computer  programs  as  functions  of  temperature.  Data  in  tabular 
form  was  read  into  the  computer  and  whenever  a  value  between 
entries  was  desired,  it  was  obtained  by  linear  interpolation.  The 
temperature  distributions  were  read  into  the  computer  as  two-dim¬ 
ensional  tables  and  linear  interpolation  was  again  used  for  inter¬ 
mediate  values.  In  all  cases,  the  tables  were  constructed  so  as 
to  tend  to  minimize  errors  arising  from  the  use  of  linear  inter¬ 
polation. 

Figure  35  is  the  comparison  between  the  experiments  and  the 
"predictions"  obtained  from  the  rracture  algorithm.  These  curves 
were  made  using  the  computer  programs  described  in  Appendices  II 
and  III  in  conjunction  with  the  material  property  curves  described 
earlier  in  this  section.  The  agreement  in  both  cases  between  the 
experiments  and  the  theory  employing  the  tension  data  is  fairly 
remarkable  in  that,  to  our  knowledge,  no  other  investigators  have 
achieved  as  good  an  agreement  between  two  sets  of  experiments  - 
much  less  three  sets  as  is  the  case  here. 
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Probability  of  Failure  i/  N-H  Probability  of  Failure  i/N+l 


Figure  35  Comparison  of  Theory  and  Experiment  for  Beam  and  Disk 

(Using  Weibull  Tensile  Strength  distribution  with  3 

m  *»  3.3,  cr,®  19,700  psi,  o  »  7,000  psi,  and  V  »  0.0982in. 

U  U  l* 
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SECTION  V 

DISCUSSION  OF  RESULTS 

1,  RELIABILITY  OF  RESULTS 

We  have  taken  the  position  in  this  investigation  that  we  have 
not  proved  the  applicability  of  our  fracture  algorithm,  but,  that 
we  have  established  its  potential  for  predicting  therraal/mechanical 
response.  To  justify  this  viewpoint  in  the  light  of  the  amazingly 
close  predictions,  we  are  obliged  to  expose  a  few  of  the  scientific 
shortcomings  of  our  program.  The  following  unanswered  questions 
all  raise  doubts  concerning  the  veracity  of  our  results. 

(a)  Have  all  the  materials  used  in  the  bean._ ,  disks, 
and  ten"ion  specimens  been  drawn  from  the  same 
statistical  population? 

(b)  Is  high  purity  aluminum  oxide  a  series  material? 

•  (c)  Temperature  dependent  elastic  and  thermal  properties 
and  statistical  strength  parameters  were  taken 
from  the  literature.  Do  they  represent  the  material 
used  in  this  investigation? 

(d)  What  is  the  magnitude  and  character  of  the  various 
parasitic  stresses  present  in  the  different  test 
specimens? 

(e)  Was  the  material  used  in  our  specimens  homogeneous; 
for  example,  is  the  surface  and  volume  material 
identical? 

Without  dwelling  on  the  desirability  of  resolving  these 
points,  let  us  say  that  a  more  ambitious  program  should  address 
itself  to  such  questions. 

2.  SENSITIVITY  OF  RESULTS 

Due  to  the  small  disparity  between  the  predicted  and  measured 
response  of  the  beam  and  disk  prototype  structures,  it  is  now 
unnecessary  to  conduct  the  error  analysis  originally  contemplated 
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at  the  beginning  of  the  program.  Instead,  a  variation  of  param¬ 
eters  study  was  performed  which  yields  more  general  information. 
Specifically,  variations  of  plus  and  minus  10  percent  (+  10%)  were 
introduced  into  each  material  property  versus  temperature  curve 
entering  into  the  response  analyses.  The  resulting  theoretical 
curves  are  compared  with  the  data  in  Figures  36  through  41. 

Referring  to  these  figures,  it  is  readily  observed  that  the 
beam  is  virtually  insensitive  to  the  mechanical  properties  e^(T), 
E(T)  and  v(T).  Poisson's  ratio  was  recognized  to  be  unimportant 
in  the  stress  analysis  of  the  beam.  However,  the  fact  that  the 
thermal  strain  and  modulus  of  elasticity  do  not  play  a  large  role 
in  the  case  of  the  beam  is  probably  due  to  the  fact  that  the 
temperature  distribution  is  too  linear.  With  respect  to  the 
Weibull  parameters  m,  ou  and  oQ  the  beam  is  observed  to  behave 
fairly  sensitively. 

The  effect  of  Poisson's  ratio  was  included  in  the  stress 
analysis  of  the  disk.  Consequently,  it  is  interesting  that  the 
+  10%  variations  in  v(T)  do  not  significantly  alter  the  theoretical 
curve  in  Figure  38.  The  other  figures  indicate  that  the  disk  is 
fairly  sensitive  to  m  and  very  sensitive  to  b^CT),  E(T),  ou  and 
o0<  It  is  also  interesting  to  note  that  the  disk  curves  for 
+  10%  6t(T)  and  +  10%  E(T),  Figures  36  and  37,  are  identical. 

This  would  be  explained  if  the  parameters  always  occurred  as  a 
product,  e^E,  as  is  the  case  in  the  stress  analysis  of  the  beam 
Equation  (23).  That  this  is  not  the  situation  in  the  case  of  the 
disk  is  readily  observed  by  inspecting  Equations  (25)  through  (34) 
for  the  stress  analysis  of  the  disk.  A  simple  example,  however, 
will  clear  things  up. 
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Probability  of  Failure  i/  N-H  Probability  of  Failure  of  i/N-H 


Figure  36  Effect  of  10  Percent  Variation  in  Thermal  Strain 
Upon  the  Predicted  Cumulative  Distributions  for 
the  Beam  and  Disk 


Probability  of  Failure  i/N+l  Probability  of  Failure  i/N  +  l 


Figure  37  Effect  of  10  Percent  Variation  in  Modulus 

of  Elasticity  Upon  the  Predicted  Cumulative 
Distributions  for  the  Beam  and  Disk 


74 


Probability  of  Failure  i/N  +  l  Probability  of  Failure  i/ N+l 


Figure  38  Effect  of  10  Percent  Variation  in  Poisson's  Ratio 
Upon  the  Predicted  Cumulative  Distributions  for 
the  Beam  and  Disk 


75 


0  1  2  3  4  5  6  7  8  9  10  II  12  13  14 

Time  toFailure(minutes) 


Figure  39  Effect  of  10  Percent  Variation  in  Weibull  Parameter 
Upon  the  Predicted  Cumulative  Distributions  for  the 
Beam  and  Disk 
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Figure  40  Effect  of  10  Percent  Variation  in  Weibull  Parameter 
Upon  the  Predicted  Cumulative  Distributions  for  the 
Beam  and  Disk 


Probability  of  Failure  i/  N-H  Probability  of  Failure  i/  N-H 


Figure  41  Effect  of  10  Percent  Variation  in  Weibull  Parameter  a 
Upon  the  Predicted  Cumulative  Distributions  for  the 
Beam  and  Disk 
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Consider  two  singly  fitting  concentric  rings  of  the  same 
material.  Assume  that  when  both  rings  are  at  the  reference 
temperature  T  ,  the  rings  are  free  from  stress.  Also,  assume 
that  both  E.J,  and  F  are  functions  of  temperature.  Now,  let  us 
examine  what  happens  when  the  inside  ring  (number  1)  is  heated 
to  and  the  outside  ring  (number  2)  is  heated  to  T2.  If  the 
thickness  of  the  ring  6  is  small,  the  radius  of  both  rings  may 
be  taken  as  R  and  thus,  after  being  heated  to  T-^  and  T2 
respectively,  the  stresses  in  the  rings  can  be  computed  from 
elementary  theory  to  be 

C  £t(Ti)  ~  eT^2^ 

1  " "  t  ~  ~ — r~ 

E(Tp  +  EXT^y 
ff2  ~  ”  al 

The  form  of  Equation  (35)  demonstrates  that  the  effect  of 
a  +  N7>  variation  in  e^(T)  will  be  identical  to  a  +  N70  variation 
in  E(T) 

3.  FUTURE  RESEARCH 

Almost  all  of  the  work  which  has  been  done  in  the  field  of 
statistical  fracture  theory  has  concerned  itself  with  the  analysis 
problem  as  contrasted  with  the  design  problem.  The  present  study 
is  no  exception  since  we  attempt  to  predict  the  behavior  of  a 
known  element  under  a  specified  environment.  If  we  are  to  deal 
effectively  with  the  problems  of  how  to  modify  an  inadequate 
structure  or  how  to  proportion  a  brittle  component  for  minimum 
weight  when  operating  under  severe  temperatures,  we  must  address 
ourselves  directly  to  the  problem  of  design.  Our  fracture  algorithm 
provides,  to  be  sure,  a  basic  tool  for  design;  but,  it  does  not 
constitute  a  design  procedure. 

To  improve  our  analysis  capabilities,  there  are  two  situations 
to  consider.  In  the  first,  we  must  extend  our  capability  so  that 
we  may  predict  the  performance  of  simple  materials  (linear, 


(35) 

(36) 
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homogeneous,  isotropic,  series)  in  more  severe  temperature 
environments.  In  the  second,  we  must  try  to  characterize  the 
behavior  of  more  complicated  materials  (nonlinear,  anisotropic, 
parallel)  under  room  temperature  and  steady  state  elevated 
temperature  environments . 
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APPENDIX  I 


.  STATISTICAL  FRACTURE  THEORY  FOR  COMBINED  STRESS  CONDITIONS 
1.  INTRODUCTION 

In  his  classic  paper  of  1939  (Reference  14),  Weibull  develop¬ 
ed  an  expression  for  the  fracture  probability  of  a  brittle  material 
under  a  polyaxial  stress  state.  Using  a  different  point  of  view, 
we  shall  expand  on  his  brief  statistical  treatment  of  this  com¬ 
bined  stress  problem  and  we  shall  extend  our  results  to  cases 
of  varying  mechanical  and  thermal  loading  and  to  materials  which 
cannot  be  represented  by  the  Weibull  distribution  function. 

Briefly,  it  is  our  objective  to  establish  a  fracture  surface, 
i.e.,  to  find  a  relationship  among  the  strengths  achieved  under 
various  stress  states.  The  usual  approach  to  this  problem  in 
either  brittle  or  ductile  materials  is  to  find  a  property  common 
to  all  stress  states  that  will  indicate  failure  or  nonfailure. 

In  ductile  materials  the  distortion  energy  represents  such  a 
property  since  incipient  flow  occurs  in  any  stress  state  in  which 
the  distortion  energy  is  equal  to  the  distortion  energy  obtained 
in  a  tension  specimen  at  yield.  Stated  in  another  way,  we  can 
correlate  yielding  under  any  stress  state  with  the  distortion 
energy.  Our  approach  for  brittle  materials  is  completely  analo¬ 
gous  -  we  shall  try  to  find  a  property  that  will  correlate  with 
the  reliabilities  associated  with  the  various  possible  combined 
stress  conditions. 

To  avoid  the  "size  effect"  problem  observed  in  the  strength 
of  brittle  elements,  i.e.,  increasing  fracture  stress  with  de¬ 
creasing  volume,  we  shall  confine  our  study  to  finite  volumes 
AV  of  fixed  size.  We  assume  that  both  the  material  and  the 
stress  state  in  these  basic  volumes  are  homogeneous  and  that  the 
materials  used  in  all  the  volumes  to  be  considered  have  been 
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drawn  from  the  same  population.  In  addition,  we  shall  re¬ 
strict  the  study  to  brittle  materials  that  are  statistically 
isotropic,  i.e.,  the  distribution  of  strengths  obtained  from 
an  indefinitely  large  number  of  basic  volumes  will  be  identi¬ 
cal  in  every  direction. 

We  shall  assume  that  the  principal  stresses  S-^,  S2, 
which  act  on  a  basic  volume  are  proportional  to  a  load  factor 
S,  i.e., 


S1  = 
S2  * 


as 

es 

7S 


(37) 


where  0,3,7  are  constants  which  define  the  stress  state.  Then, 
the  strength  of  a  basic  element  will  be  taken  as  the  maximum 
load  factor  that  it  can  equilibrate.  Failure  of  the  element 
is  represented  by  its  inability  to  equilibrate  the  applied 
loading.  It  is  important  to  point  out  that  it  is  possible 
for  cracks  to  initiate  and  propagate  within  the  basic  volume 
without  causing  its  failure.  Materials  in  which  cracks  can 
be  arrested  or  which  provide  alternative  load  paths  when 
local  failures  occur  are  classified  as  parallel  or  series- 
parallel  materials.  If  a  local  failure  necessarily  leads  to 
overall  failure,  the  associated  material  is  called  a  series 
or  "weakest  link"  material.  One  can  advantageously  adopt  an 
infinitesimal  basic  volume  for  the  series  material  and,  as 
we  shall  subsequently  discuss,  combined  stress  testing  is 
greatly  simplified  in  this  case. 

Only  the  tensile  or  cohesive  mode  of  failure  will  be 
considered  in  this  investigation.  We  shall  ignore  the  in¬ 
fluence  of  compressive  or  shear  stresses  on  the  strength  of 
a  brittle  material.  The  potential  usefulness  of  this  tension 
criterion  is  a  consequence  of  two  observations;  first,  that 
the  shear  strength  of  brittle  materials  is  usually  an  order 
of  magnitude  greater  than  the  tensile  strength  and,  second, 
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that  it  is  extremely  di_ficult  to  eliminate  tensile  stresses 
from  prototype  or  laboratory  elements.  Almost  every  structural 
failure  of  a  brittle  component  can  Ke  attributed  to  the  pres¬ 
ence  of  some  distribution  of  tensile  stresses. 


2.  TWO-DIMENSIONAL  THEORY 


Heuristic  Development-Single  Loading 


When  we  attempt  to  describe  the  statistical  fracture 
strength  of  a  finite  volume  of  material  under  a  uniaxial  stress 
state,  the  axial  stress  (strain)  is  the  only  reasonable  choice 
for  the  statistical  variate.  Taking  a  general  form  for  any 
cumulative  distribution  function,  we  can  write  the  fracture 
probability  F  for  the  uniaxial  stress  state  as 


F(s)  -  1  -  exp  [-  g(o)j 


where  AV  is  the  specified  volume  of  the  basic  element,  v  is 
a  volume  of  unity  and  a  is  the  axial  stress.  The  delineation 
of  the  constant  AV/v  does  not  affect  the  generality  of  this 
expression  and  in  the  special  case  of  a  series  material  it  pro¬ 
vides  a  convenient  representation.  If  we  examine  the  strength 
of  a  basic  volume  of  an  isotropic  material  under  a  general  homo¬ 
geneous  stress  state,  it  follows  that  failure  will  depend  only 
on  the  three  principal  stresses  acting  on  the  unit.  Thus,  the 
probability  of  failure  of  the  basic  volume  can  be  designated 
as  F(S^,  $2,  S^)  where  the  three  principal  stresses  are  taken 
as  the  statistical  variates.  For  this  case  we  shall  take 
Equation  (38)  in  the  form 


■log  fl-F(SlfS2lS3)’ 


0  g(s1,s2,s3) 


For  a  specified  reliability  (1-F),  Equation  (39)  becomes  gCS^.S^.S 
equals  constant,  which  defines  our  fracture  surface. 

On  the  basis  that  failure  is  caused  only  by  tensile  stresses, 
it  seems  reasonable  to  look  for  the  function  g  within  the 
collection  of  all  possible  tensile  stresses  which  can  occur  at 
any  point  in  the  basic  voLurae.  In  the  plane  stress  problem,  we 
cau  relate  the  normal  stress  crn  acting  in  any  direction  to  the 
principal  stresses  through  the  expression 
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(40) 


2  2 

a  =  S-,  cos  0  +  S0  sin  9 
a  1  L 

where  9  is  the  angle  between  an  and  S-^.  As  9  sweeps  through 
all  values  from  -  r /2  to  r  /2,  Equation  (40)  describes 
every  possible  normal  stress  acting  at  a  point.  The  normal 
stresses  associated  with  the  various  directions  described  by 
9  are  shown  in  Figure  42  for  several  different  stress  states. 

The  question,  now,  is  what  are  the  distinguishing  features  of 
these  figures  which  will  reflect  the  differences  they  cause  in 
a  material's  response? 

The  most  obvious  first  guess  is  to  differentiate  among 
these  stress  states  by  comparing  the  areas  associated  with  the 
tensile  normal  stresses.  However,  this  approach  does  not  re¬ 
flect  the  possibility  that  the  magnitude  of  the  stresses  may 
have  a  different  influence  than  their  extent  or  distribution. 

For  example,  hydrostatic  tension  and  pure  tension  stress  states 
which  lead  to  the  same  area  are  depicted  in  Figure  42a  where 
we  observe  that  one  peak  stress  ia  twice  the  other.  Experience 
indicates  that  the  pure  tension  state  is  the  moi’t  critical.  On 
the  other  hand,  when  the  maximum  hydrostatic  tension  and  pure 
tension  stresses  are  equal,  as  indicated  in  Figure  ^3b  ,  our 
intuition  would  select  the  hydrostatic  state  as  the  more  critical. 
Unlike  pure  tension,  oriented  flaws  such  as  cracks  cannot  avoid 
exposure  to  high  normal  stresses  by  assuming  a  preferred  direc¬ 
tion  since  all  directions  experience  the  same  stress  under 
hydrostatic  tension.  This  implies  that  a  maximum  stress  theory 
is  inapplicable  and  unc.onservative,  and  indeed,  evidence  exists 
to  support  this  contention  (Reference  15) . 

The  two  examples  depicted  in  Figure  43  suggest  that  we 
"weight"  the  ordinates  of  the  normal  stress -theta  diagrams  for 
different  stress  states  and  then  compare  their  areas.  Assuming 
a  statistically  isotropic  material,  the  weighting  should  be  in¬ 
dependent  of  the  orientation  of  the  normal  stress,  0.  We  might 
use,  for  example,  a  power  function  to  modify  the  normal  stresses, 
i.e. , 
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■Pure  Tension; Of,= Of,  (2S,0) 


Hydrostatic  Tension 
Op,=CTn(S,S)  r 


-7T/2  -7T/4  7T/4  TT/2  9 

(a)  Equal  Areas  -Unequal  Maximum  Stresses 


Pure  Tens  ion :  On5  Of,  (S,0)- 


Hydrostatic  Tension  :Of,=Of,(S,S) 


"77/2  TT/2  & 

(b)  Equal  Maximum  Stresses -IJnequa1  Areas 


Figure  43  Comparison  of  Hydrostatic  Tension  and  Pure 
Tension  States 
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(41) 


D  a' 


n 


where  D  and  k  are  constants.  This  alteration  results  in 
the  dashed  curve  shown  in  the  left  side  of  Figure  44  .  if  the 

normal  stress  distribution  for  several  stress  states  are  weighted 
in  this  fashion,  we  could  compare  the  areas  of  the  resulting 


curves,  i.e.. 


g(S-^ ,  S2)  =  Area  =  D 


°n  d9 


(42) 


0  >0 
n 


where  the  integration  extends  over  those  values  of  9  where 
the  normal  stress  is  non-negative.  Because  of  symmetry  we  need 
consider  only  the  positive  normal  stresses  in  the  interval  zero 
to  tt/2.  To  account  for  the  possibility  that  tensile  stresses 
below  a  certain  magnitude  a ^  may  not  cause  failure,  we  may 
choose  to  weight  the  difference  (<?n“Op  as  shown  in  the  right 
half  of  Figure  4i.  The  associated  area  is  given  by 


g(s1}s2) 


=  Area  =  D 


I 


(a. 


n 


k 

Cg)  d9 


(43) 


V°Jt 

Certainly,  the  use  of  a  power  function  to  weight  the  normal 
stress-theta  diagrams  is  completely  arbitrary  and  there  are 
many  other  ways  of  manipulating  and  distorting  such  curves.  Our 
problem  is  to  find  a  weighting  function  that  will  reflect  the 
influence  of  stress  state  on  the  reliability  of  a  basic  volume. 
Denoting  the  weighting  function  by  f,  the  fracture  probability 
becomes 

tt/2 

-  AV 


F(S^,S2)  3  1  -  exp 


1 


£(an) 


dQ  (44) 


We  are  now  in  a  position  to  describe  certain  guidelines 
for  the  selection  of  f.  First,  to  account  for  the  possible 
existence  of  a  zero  fracture  probability  stress  a-, ,  we  must  take 


aa 


f  =  f<an  ‘  ai> 


an- 


(45) 


0 

f  =  0  ans  at 

The  latter  condition  implies  that  both  S^s,  and  S2~  anc^ 
that  in  such  cases  F  =0.  At  the  other  extreme,  we  expect  that 
fracture  is  a  certainty  when  either  or  S2  is  positive  and 
unbounded;  hence,  F  =  1  implies  that 

f  -*•  00  when  S^,S2  +  « 

Furthermore,  we  would  expect  on  physical  g...  mds  that  the 
failure  probability  would  increase  continuously  with  increasing 
principal  tensile  stresses;  thus, 

f. . .continuous  and  monotone  increasing. 

Finally,  f  must  be  chosen  in  such  a  way  that  the  associated 
FCS-^.S^)  fits  the  cumulative  distribution  curves  obtained  from 
fracture  tests  conducted  using  various  stress  states.  In  parti¬ 
cular,  it  is  necessary  that  fracture  data  obtained  under  pure 
tension  be  represented  by  F(S^,0)  or  F(0,S2)  and  that  hydro¬ 
static  tension  data  be  represented  by  F(S^,S^),  This  is  a 
standard  problem  in  curve  fitting  and  one  proceeds  by  selecting 
a  reasonable  and  versatile  form  for  f  which  contains  0*  and 
n  additional  parameters  a^;  i.e.,  £  »  f  [(cr^  -  a ^);  a^,a2,...anj 
These  parameters  are  chosen  so  that  the  curve  for  F  passes 
"as  close  as  possible"  to  each  data  point.  For  a  aeries  material 
we  note  that  the  parameters  which  provide  an  exact  fit  to  an 
infinite  amount  of  data  are  intrinsic  phenomenological  strength 
properties  of  the  material.  Otherwise,  they  characterize  ap¬ 
proximately  the  strength  of  the  basic  volume. 

The  following  typical  functions  may  be  useful  candidates 
for  f(S^, 82,22): 


For  a  Co, : 
n  i 


f  =  0 


For  a  2a. : 
n  1 


a  -  a„  k 


a  -  a , 


(-^-)  ■  ;  k>  0  ; 

c 

**  c-V-^ 

c 

(46) 

Z  +  a1>0  ; 

k  ^  0 

(4  7) 

9*  slqZ  ■+■  "t  ... 

+  a„Xn  ;  ai2  0 

(48) 

exp  1-1 

(49) 

exp  1k  -  1  ;  k  >  0 

(50) 

exp  (exp  1  -  1)  -  1 

(51) 

where  k,  a^,  ac,  and  the  a^'s  are  statistical  parameters  and 

where  or  is  the  normal  stress.  The  first  of  these  functions, 
*-k  n 

f  "  X  ,  is  the  one  adopted  by  Weibull  (Reference  14) .  We  shall 
interpret  f  as  the  generalized  normal  stress  and  T  as  the 
generalized  normal  stress  vector;  in  two  dimensions  7  «  7(£,9), 

h.  Specific  Formulas 

In  two  dimensions  the  fracture  probability  is  given  by 

F(S1,S2)  o  1  -  exp  -  j  f(on  -  o,)  d©  !  (52) 

o 

where  the  integration  extends  over  the  region  where  0 
specifically. 


7  2 

o,  84  S-.  cos“9  +  S~  sin  9>o, 
n  L  l  £ 


(5J) 


Theta  is  the  angle  between  a  and  S, , 

11  i> 
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This  inequality  leads  to  three  distinct  cases: 

(1)  S1>S2>ct,>0: 

So  "  ^2 

(2)  S^~  a* ;  S2  -  a ^ : 


0Q  =  cos 


I  /oT  -  S0 

,1\  x  2 


S1  "  S2 


(3)  SL<  Ojj ;  S2  < Og 


9q  =  0  (F  -  0) 

Selecting  the  Weibull  form  for  f  given  by  Equation  (46) 
we  shall  determine  the  failure  probability  for  pure  tension  and 
for  hydrostatic  tension. 

Pure  Tension:  =>  ;  S2  =  0 


FCcO  »  1  -  exp 


cos  " 1  W°t  n 

aii  I  atcos^9  ~  o;  k  ! 

~  (-^—z - -)  ue:  (54) 


When  *  0,  F(at)  becomes 

f  av  Vtt  F(k  +  -4“) 

F(ot)  «  1  -  ekp  I  -  —  ~T~\mr  (~)  ( 

*  c 


55) 


Hydrostatic  Tension:  “  S2  “  >o 

F(bh)  -  1  -  exp 


hV  Jt  ,  % _ % 

r  t  1  o..  } 


k 


(56) 
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This  expression  has  the  same  form  as  the  Weibull  distribution 
function;  consequently,  the  estimation  of  aQ,  and  k  can 

follow  well  developed  procedures  for  this  purpose  (References  2, 
14,  and  16).  If  biaxial  hydrostatic  tension  data  can  be  obtained 
the  parameter  determination  for  all  of  the  forms  indicated  in 
Equation  (46)  through  (51)  will  be  greatly  simplified  since  all 
the  f's  become  independent  of  theta. 

c.  Multiple  Loadings 

If  our  basic  volume  is  subjected  to  a  number  of  different 
loadings,  we  shall  attempt  to  correlate  fracture  with  the 
"worst"  conditions  that  can  arise.  It  is  assumed  that  stresses 
which  act  on  the  basic  volume  prior  to  fracture  do  not  signi¬ 
ficantly  effect  the  resistance  of  the  material.  Now  then,  we 
shall  separately  consider  every  possible  direction  in  our  basic 
volume  and  identify  the  largest  normal  tensile  stress  which 
acts  throughout  the  load  history;  i.e.,  &n(t)  .  Eacn  of  these 

maximum  normal  stresses  will  be  weighted  to  form  the  maximum  gen¬ 
eralized  normal  stress  fmax  which  will  then  be  summed  over  theca 
in  the  usual  way  to  form  g.  For  example,  consider  three  distinct 
loadings  which  give  rise  to  the  same  principal  directions  such 
as  illustrated  in  Figure  45a.  The  collection  of  maximum  normal 
stresses  is  indicated  as  the  envelop  of  the  three  solid  curves 
and  the  weighted  envelop  curve  is  shown  as  a  dashed  curve.  We 
note  that  symmetry  is  preserved  across  the  on  and  the  9  axes.  The 
fracture  probability  may  then  be  written  as 


where  we  recall  that  the  constants  o,  and  the  a.'s  are  statisti 

&  L 

cal  parameters. 

Although  we  have  chosen  in  our  previous  work  to  represent 
the  o  -  9  relationship  in  cartesian  coordinates,  there  is  no 
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fundamental  reason  for  doing  so.  For  the  multiple  loading  prob- 
.lem  it  is  more  convenient  to  use  a  polar  representation  as 
shown  in  Figure  45b  where  the  normal  stress  vector  an  and  the 
generalized  normal  stress  vector  7  are  represented  in  a  natural 
way.  The  x-y  axes  are  assumed  to  be  fixed  in  the  body  and  both 
the  principal  stresses  and  directions  are  taken  as  time  dependent, 
Calculating  the  area  under  the  weighted  maximum  normal  stress 
envelop,  fmaxj  to  establish  g,  the  fracture  probability  in  two- 
dimensional  polar  coordinates  becomes 

F  '  °xx(t)>  =  1 


_  .£ 


yy 


XV' 


or 


exp 


r 


r  ft 

\  j-AV  | 


f  max  ^ 


J 


r  dr  d9 


o  o 


F  [®xxW>  °yy(c)*  V;cb  =  1-**P  dr  J  4*  £2  { j“x  ®n(t)J 


-  n 


•  ai}  d0) 


(53) 


where  we  note  that  symmetry  is  available  only  through  the  origin 
and  where  the  normal  stress  is  given  by 

**  Vx(c)  cos~^  +  t)  sin"9  +  Oj-yCfc)  Sin29  (59) 


Comparing  Equations (53)  and  (57), we  find  that  the  polar  re- 

7 

presentation  integrates  f“  rather  than  f  over  theca.  This  does 
not  lead  to  an  essential  difference  in  the  fracture  probability 
expressions  since  in  each  case  we  muse  find  an  f  that  fits  our 
data.  I’or  example,  if  f  is  given  by  Equacion  (4b)  the  ex¬ 
ponent  is  absorbed  by  the  parameter  k  and  differences  in  multi¬ 
plicative  constants  are  absorbed  by  <r  . 

d.  Thermal  Loading 

It  has  tacitly  been  assumed  throughout  this  study  that  under 
fixed  environmental  conditions  (temperature,  atmosphere,  humidity) 
the  cumulative  strength  distributions  for  our  basic  volume  do 


not  vary  with  time.  We  do  not  account,  for  example,  for  the 
possible  effects  of  static  fatigue,  creep,  or  continuous 
chemical  changes  in  our  material.  Where  this  assumption  holds 
we  can  associate  a  unique  function  f  with  every  environmental 
condition.  Specifically,  for  a  uniformly  distributed  temperature 
T  in  our  basic  volume,  we  shall  designate  the  generalized  normal 
stress  by  £p(<?n*aiJ  >  a^)  •  In  practical  situations  one  usually 
tries  to  use  as  few  different  forms  of  f  as  possible.  For  a 
given  form,  however,  the  statistical  parameters  must  be  taken 


temperature 

dependent. 

Thus,  fT, 

can  be  written, 

fT  -  q 

k- 

Ojtt) 

;  »tCT)] 

TqS  T  Tl 

=  f? 

k- 

a2(T) 

;  bj(T>] 

h-  T“T2 

(60) 

f ®  f 

T  m 

l°n  ' 

°m<T> 

* 

I  T*T 

m- 1-  *•  m 

where  the  temperature  range  of  interest  is  T  ^  T*;,  T  and  where 
the  symbols  Og,  om,  b^,  c^.  represent  temperature 

dependent  statistical  parameters. 

When  a  basic  volume  is  subjected  to  a  temperature  history 
T(t)  and  a  stress  history  o^(t),  we  once  again  identify  the 
"worst"  condition  in  every  direction.  Here,  we  should  recognize 
that  the  largest  normal  tensile  stress  in  a  given  direction  does 
not  necessarily  lead  to  the  largest  generalized  normal  stress. 
What  we  must  find  is  the  maximum  combination  of  normal  stress 
and  weighting.  For  a  single  stress  state,  the  envelop  of  maxi¬ 
mum  generalized  normal  stresses  is  given  by 

f  “  max  £  [  cr  -  o.(T);  a.(t)l  (6L) 

max  *  n  i.  i  i 

When  the  stress  state  varies,  this  maximum  envelop  is  defined  by 
fma*  "  in?x  fT(c)  *  °n(t)-  »  Jil  (  o2) 
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where  the  statistical  parameters  and  a^  depend  explicitly 
on  the  temperature  and  implicitly  on  the  time  parameter  t. 
Adopting  polar  coordinates  with  fixed  axes  x-y  in  the  basic 
volume,  we  form  the  function  g  by  summing  this  f  .  over 
theta;  thus, 

F  [°xx(c)>  V(t)l  '  l‘exp 

o 

*1 

•V2 

/on(C>  -  <j^[t(C)  ;  at  ;'t(c)'V  (63> 

’  i 

where  f^,  is  defined  by  Equation  (60)  and  on  is  given  by 
Equation  (59).  This  expression  predicts  the  fracture  proba¬ 
bility  of  a  basic  volume  under  biaxial  time  varying  stresses 
and  subjected  to  a  changing  temperature  environment. 


7T 


<  - 


AV 

v 


1  /. 


\mtX  £T(t) 


3.  THREE-DIMENSIONAL  THEORY 


a.  Single  Loading 

Following  the  two-dimensional  development,  we  propose  to 
correlate  the  behavior  of  the  basic  volume  under  a  general 
stress  state  with  some  function  of  the  totality  of  normal  ten¬ 
sile  stresses  in  three  dimensions.  In  the  belief  that  any 
tensile  normal  stress  can  occasionally  cause  a  fracture  by  vir¬ 
tue  of  its  magnitude  and  relative  orientation  with  respect  to 
a  random  flaw,  we  once  again  assume  isotropy  and  weight  each 
normal  stress  vector  on  according  to  its  magnitude  only.  Using 
the  polar  coordinates  defined  in  Figure  46a,  the  normal  stress 
in  three  dimensions  can  be  written  in  terms  of  the  three  princi¬ 
pal  stresses;  thus, 

oa  =  cos4^  (S^  cos~tf  +  S,  sin~tf  )  +  sin^d  (64) 

The  normal  stress  diagrams  associated  with  several  important 
stress  states  are  sketched  in  Figure  46  where  we  observe  the 
following: 

(1)  Symmetry  makes  it  possible  to  consider  only 
one  quadrant, 

(2)  The  hydrostatic  stress  state  gives  rise  to  a 
spherical  surface. 

(3)  The  zero  probability  stress  a*  can  be  represented 
as  a  sphere  of  radius  o,  which  can  be  subtracted 

A 

from  each  diagram  to  yield  dn  -  a.,  . 

(4)  If  the  three  principal  stresses  are  positive, 
tne  volume  of  the  normal  stress  diagram  is 
given  by 

V  “  Tpjj  1 5(S  jj)  +  3(S£S.?+S£j2‘*‘S^S^*‘S2SybS2S|^s-S^S-})  +  2S^S-,S-jj 

(65) 
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Figure  46  Normal  Stress  Diagrams,  Three-Dimens ional 


This  expression  is  useful  as  a  computer  check. 


Now,  adopting  the  same  types  of  weighting  functions  f 
described  in  Section  2-a,  we  can  write  an  expression  for  the 
fracture  probability  of  a  basic  volume  subjected  to  a  three- 
dimensional  stress  state;  thus, 


F(S-j_,S2>S2)  =  1  -  exp  { - ^  JJj  rZ  dr  cos«5  df\  (66) 


where  the  integration  extends  over  that  portion  of  the  first 
quadrant  where  the  normal  stress  is  greater  than  the  zero 
probability  stress.  More  specifically, 


tt/2 


F(S1,S2,S3)»  1-exp  d *J  -r  f3  [vCTr  aJ 


costf  d^ 


where  we  delineate  four  distinct  cases: 


(1)  SL  ;  S2  tS3 


*1 

d. 


o 

7l/2 


(2)  8^82*83^ 

<6t  “  0  1 

L  >  F  13  0 

^  »  oj 

(3)  S^Sj-Oj,  S^a- 


«j  -  s3 


cos  1 - — - — — ->r 

S^cos  v  +  S9sin“^- 


- 


(67) 
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(4)  S3?^  ,  S ^  S  2  £  ^ 

=  cos  1 


S3  gl 

SyS-^cos2^  -  S2sin^ 


JU 


=  7t/2 


The  integration  limits  in  cases  3  and  4  are  derived  in  Figure  47. 

Using  the  Weibull  form  for  f  given  by  Equation  (46),v7e  shall 
determine  the  failure  probabilities  for  several  important  cases. 

Pure  Tension: 

Case  4;  <?n  =  ^(0,0,^);  </>L  =  cos"1  -  vr/2 


?r/2  7t/2 


F(0,0,crt)  =  1-exp  (  -  ^  /  d*  /  ± 


!  ■  ^j-sin2^-^  j 


3k 


I 


cos«$  d«5  )  (63) 


0  cos’1  '/  (cr^-c-jg ) /crt 


For  <rtf  «  0, 


g(0,0,Pt)  -  6(6k+l)  (  or* 


3k 


(69) 


Comment:  This  formula  for  g  is  much  simpler  than  the  correspond¬ 
ing  g  in  the  two-dimensional  case  described  by  Equation  (55).  It 
has  the  same  form  as  that  found  in  the  conventional  Weibull  distri¬ 
bution  function  and,  consequently,  the  estimation  techniques  de¬ 
scribed  in  References  2,  14  and  16  can  be  used  to  establish  the 
statistical  parameters  cr  and  k. 

V 
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Hydrostatic  Tension:  3-D 


Case  1;  \  ?  ip  ;  ^  —  0  ;  -  71  / 2 


U 


f 


77  /2  Tr/2 


r 


F^sh’sh’V  = 


l-exp  /  -  —  j  drt 


,S  ~a  '■  3k  i 

1  '  sh  ai?  v 

i  i  ~ — “I  cosrfdd) 

i  'c  i  ! 


0  0 


(70) 


-  ;  ®h_cr7  • 
g(s^sh>sh)  -  6 

V  c  / 


3k 


(71) 


Comment:  We  can  compare  the  magnitudes  of  hydrostatic  tension 
and  pure  tension  when  each  produces  the  same  failure  proba¬ 

bility  in  a  basic  volume  AV.  Equating  the  two  functions  g  given 

* 

by  Equations  (69)  and  (71)  and  taking  =  0,  we  obtain 


Is  h]  IjJiU*. 

6k+l  j 

s  l 


rr 

\  t 


(72) 


Hydrostatic  Tension:  2-D 

Case  3;  on  =  >  °i  "  0  5  3  0  '  =  1T^2 

tt/2  rr/2 


F(°rh,CTh,°^  “  l~Q*? 


L 


r  '  2 

/  ,  /a.  cos  «J1 

dtf  j  ^  - -  cos«J  dd 

j  \  c  l 

0 


(73) 


nx  rr3/2  V(3k+l)l\\ 

h  h  12(3k+£)  r(3k+j)  \ffc/ 


3k 


(74) 
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Comment:  The  expression  for  g((Th,ch,0)  is  more  complicated  than 
the  corresponding  g(o,^,o^)  obtained  from  Equation  (56)  in  our 
two-dimensional  treatment  of  the  hydrostatic  tension  plane  stress 
problem.  For  equal  reliabilities  we  can  compare  the  biaxial  hydro¬ 
static  state  with  the  pure  tension  state  by  equating  the  functions 
g  given  by  Equations  (69)  and  (74);  thus, 


\ 


r(3k+i)  1 
r"  (3k+i)  | 


l/3k 


3-D  Theory 


(75) 


We  form  the  corresponding  expression,  using  the  two-dimensional 
theory,  by  equating  Equations  (55)  and  (56)  and  replacing  the 
letter  k  by  k' ;  thus, 


>D  Theory 


(76) 


To  compare  the  (2-D)  and  (3-D)  theories,  the  same  input  data  must 
be  used  to  establish  the  respective  parameters.  Therefore,  if  each 
theory  is  fit  to  the  same  tension  data,  we  can  equate  Equations  (55) 
and  (68)  to  obtain  the  relationships  among  the  various  parameters, 


i.e 


*  J 


y;r(k'+i) 

[V 

k'  . 

7T  /  °t  \ 

3k 

2  r(k’+i)  | 

14 

6(6k+i) 

(77)  ) 

t 


The  equivalence  is  established  if 


k' 


3k 


a  0 


3(2U'+n  r(k4>! 

— 7;  rrnrrTj 


1/k' 


(78) 

(79) 


Since  only  one  parameter  appears  in  Equation  (75)  we  simply  replace 
3k  by  k'  to  see  that  Equations  (75)  and  (76)  are  identical.  In 
general,  the  (2-D)  and  (3-D)  theories  both  predict  the  same  fracture 
curve  for  the  plane  stress  problem  when  the  Weibull  f  is  used. 
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S1-S2-S,S3— S: 


Case  3;  cr  =  S  cos2^  ;  ",  =  0  ;  </T  =  0  ;  ^IT  =  rr/4 


r 


tt/2  tt/4 

f  f 


F  (S  ,S ,  -S)  =  1-exp  '  -  —  /  dii  \  'S-?9§M 

v  I  j  3  ;  -rc  , 

1)  "'o 


3k 


cos^ 


(80) 


k“l: 


k«2:. 


,,  I  o  \py  '  o  3 

8(s.s,-s)  ;^_2;(§-) 


g(s  ,s ,  -S)  =  j- 


'512 


6  3003  i 


\T?]  S  6 
2  •  (”) 
'c 


(81) 

(82) 


Comment:  In  situations  where  an  axial  prestress  is  brought  onto 
a  beam  or  column  member,  it  may  be  important  to  know  the  increased 
tensile  resistance  in  the  transverse  directions.  The  transverse 
integrity  without  prestressing  can  be  estimated  by  the  biaxial 
hydrostatic  tension  case.  Comparing  Equation  (74)  with  Equations  (81) 
and  (82)  we  find, 


for  k^l, 


3/2  r'(3+l){%\3lj  n  1 8VI)  ,Sv 

12(34)  ?04)\ac/  ^  I  ac 


(83) 


S  »  1.1222  cr. 


(84) 


for  U*2, 


3/2 


rw+ufV;6.  n  ( 5 12 V~:; )  {s_.' 

6  \T55T  i  V 


12(6-4)  2(6-4)  '/'c/ 


S  -  1.0594 


h 


(8b) 

(86) 


Thus,  our  axial  prestress  increases  the  transverse  strength  by 
12.3  percent  when  k“l  and  by  5.9  percent  when  k=2 . 
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S1=S2=-S , S3=S : 


Case  4:  a n  =  -S  cos2d  ;  0  ;  ^  =  tt/4  ;  ^  =  tt/2 


[  vr/2  tt/2 

F(-S,-S,S)  =  1-exp  j-  ~  [dir  j  \  l  ^os2^  )  cos ^ 


J 


d«5  I  (87) 


0  '  tt/A 


k=>l: 


g(-S,-S,S)  = 


"  is  V2-9  Is  x3 

s  l- 35— i(a-> 


(88) 


k=2 : 


g(-S,-S,S) 


1  >'835-512  VI)  ,S  x6 

6  im-i  I  J 


1 


3003  / V 

*  r 


(89) 


Comment:  It  is  possible  to  increase  the  axial  tensile  strength  of 
a  member  by  applying  a  transverse  compression  loading.  The  problem 
arises,  for  example,  in  certain  collet  type  grips  for  tension  mem¬ 
bers.  Such  grips  produce  a  radial  compression  state  which  is  pro¬ 
portional  to  the  tensile  loading.  The  resulting  increase  in  re¬ 
sistance  at  any  reliability  level  can  be  determined  by  comparing 
Equation  (69)  with  Equations  (83)  and  (89);  hence, 


for  k«l, 


6(6+1) 


(90) 


for  k-2. 


S  «  1.2924  rr 

t 


srd+iy 


/  835-512  VTl 
\  3oo3  / 


(91) 


(92) 


S  -  1.1299  crc 


(93) 


Hence,  an  increase  of  29,2  percent  in  the  tensile  strength  is  ob¬ 
tained  when  k«l;  13  percent  when  k=»2. 
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b.  Mechanical  and  Thermal  Load  History 


Following  the  development  of  Section  2-d,  we  shall  once 
again  consider  a  basic  volume  subjected  to  a  temperature  history 
T(t)  and  a  stress  history  'm^(t).  The  extension  of  our  previous 
treatment  of  this  problem  to  throe  dimensions  requires  that  we  use 
the  three-dimensional  form  for  the  normal  stress  and  the  appropri¬ 
ate  summation  of  the  maximum  generalized  stress  f  „  given  by 
Equation  (62), 

Using  cartesian  tensor  notation,  the  components  of  the  stress 
vector  T  acting  on  a  plane  with  unit  normal  n  can  be  written  in 
terms  of  the  stresses. 

Ti  "  Ji"j  ’  l'2>3  <*>> 

The  magnitude  of  the  normal  stress  vector,  "n,  is  given  by  the 
scalar  product  of  the  stress  vector  and  the  unit  normal  vector;  thus, 


°u  “  Tt  n i  ’  nj  ni 


(95) 


Interpreting  the  components  of  the  unit  normal  vector  n^  as 
the  direction  cosines  of  n_,  we  can  relate  the  n^'s  to  the  polar 
coordinates  shown  in  Figure  ^6a. 


n^  =»  cosi  cosy 

n2  “  cos^  sintf  (96) 

n-j  a  sin^ 


Substituting  Equations  (96)  into  (95)  we  obtain 
cf  ®  cos“^(cty  cos  y  +  sin“y)  +  0^3  s^n 

7  t 

+  cos  d  siu2y  +  ^73  sin2^  sinv  +  0 


31 


sin2^  cosy 


(97) 
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For  the  most  general  loading,  we  can  expect  only  that  the 
normal  stress  vectors  will  be  symmetic  with  respect  to  the  origin; 
consequently,  we  must  consider  all  of  the  normal  stresses  above 
the  1-2  plane.  On  this  basis  we  can  write  down  the  fracture  prob¬ 
ability;  thus. 


cos <i>  d«j  > 

i 

\ 


(98) 


where  on  is  given  by  Equation  (97)  and  by  Equation  (60). 

The  function  f  contained  in  the  definition  of  fj,  should  meet 
the  conditions  described  in  Section  2-a. 


4.  THE  FRACTURE  SURFACE 
a.  Exact  Theory 

For  a  specified  reliability  (1-F)  and  a  given  basic  volume 
4V,  Equation  (39)  defines  the  fracture  surface  which  describes  the 
resistance  of  the  basic  volume  under  various  stress  states.  We 
customarily  normalize  such  surfaces  by  relating  all  behavior  to 
that  under  pure  tension.  Then,  the  fracture  surface  becomes 


8(sr  S2>  S3^  _  , 

IUHTT*  "  1 


(99) 


All  points  on  this  fracture  surface  have  the  same  fracture  prob¬ 
ability;  points  falling  inside  have  lower  fracture  probabilities 
and  points  falling  outside  have  higher  ones. 
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In  the  two-dimensional  problem  the  fracture  surface  can  be 
represented  as  a  curve  defined  by 


g(sr  s2) 

-5T^r= 1 


(100) 


Specializing  this  formula  to  the  Weibull  power  function  with  «  0 
we  obtain, 


j 

[iij  cos2e  +  ! 

!sA  2 

— ;  sin20 

1 

L\"t;  1 

it/ 

d6 


'„2° 


V7  r<k4) 

2  I’(k+1) 


(101) 


We  have  already  evaluated  this  expression  for  the  hydrostatic  ten¬ 
sion  case  in  Equation  (76).  Numerical  results  for  this  case  are 
tabulated  in  Table  VI  together  with  two  other  stress  states.  We 
should  point  out  once  again  that,  when  the  power  form  of  f  is 
used,  the  stress  ratios  can  also  be  computed  from  Equation  (99) 
when  is  set  equal  to  zero. 

Table  VI 

STRENGTH  RATIOS  FOR  TOO- DIMENSIONAL  STRESS  STATES 


Stress  State 

Strength 

Ratio 

k=l  ' 

k°2 

;  U»3  ik^ 

1  i 

-■  -  - — —  -  - 

Pure  Tension  (<^,0) 

rr  /  r*  =» 

t'  t 

■  i  ; 

1 

1  1 

1  1 

[Hydrostatic  Tension  (%»ch) 

"i/'t  ■ 

;  1/2  ; 

(3/8) 1/2 

(5/16) 1/3 

1 

1 

I 

l 

Biaxial  Tension  1 

v*t  ■ 

2/3 

(12/19) 1/2 

(40/63) 1/3 

!  2 

» 

jPure  Shear  ( <*  , ,  - ? , ) 

j  b  S 

f  t?  ca 

tt/2 

<3/2)1/2 

(15rr/ 32) 1/3 

l 

!  1 

The  fracture  curve  associated  with  Table  VI  is  shown  in  Figure  43. 
We  observe  that  the  maximum  stress  theory,  max(S^,S^)  <_  coincides 
with  the  case  k=«>.  It  was  pointed  out  by  Weibull  (Ref .  14)  that 
corresponds  to  a  classical  deterministic  material  with  an  ultimate 
strength  equal  to  *  .  We  also  find  in  Figure  48  that  the  tensile 
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strength  increases  as  the  transverse  compressive  stress  becomes 
greater.  Physically,  the  compressive  stress  can  be  attributed 
with  closing  up  some  of  the  cracks  which  might  otherwise  be- 
critical  under  the  tensile  loading. 

As  the  compressive  stress  increases  indefinitely,  our  theory 
indicates  that  the  tensile  strength  becomes  unbounded.  This  is 
clearly  contrary  to  reality;  but,  it  is  not  a  surprising  result 
since  no  provision  has  been  made  in  our  theory  to  account  for 
compressive- shear  failures.  The  emergence  of  another  failure  mode 
provides  a  limitation  on  the  range  of  applicability  of  our  theory. 
The  ragged  line  in  Figure ^8  is  meant  to  indicate  such  a  limitation. 

b.  Approximate  Theory 

If  the  behavior  of' a  basic  volume  under  a  general  stress 
state  was  equivalent  to  the  behavior  of  three  basic  volumes  each 
under  a  distinct  principal  stress,  the  survival  of  the  basic  vol¬ 
ume  would  require  the  simultaneous  survival  of  each  unit.  Under 
these  conditions  the  reliability  of  the  basic  volume  would  be 
given  by 


1-F 

Q 


fl-F(S1)]  (l-F(S2)]  (i-F<S3)| 


(102) 


where  F(S)  is  the  fracture  probability  of  a  basic  volume  under 
a  pure  tensile  stress  S. 

For  a  Weibull  material  subjected  to  pure  tension,  the  reli¬ 
ability  is  expressed  as 


1-F(’t)  "  exp 


* 

r  ft  „  .■*.  \  m 

-  jl 

ft  u  1 

V 

■1 ,  I 

0  / 

L 

c  u 


't  ~  u 


(1.03) 


where  cr^,  r»  ,  m  are  the  Weibull  statistical  parameters 
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The  associated  approximate  combined  stress  theory 


( 


[l-Fa(S1,S2,S3)]  =  exp{- 


(Rr"u 


£iV  _ 

v  L l  ao  | 


m 


\m 


R,-<r  \  /R--~ 

+  |  _2_u  +  _3 — u 

cr  ,  c  i  i 

0  /  \  0  /  J 


(1( 


where 


Ri  =  Si 


when 


S.  >  o 
l  u 


Ri  a  "u  when  <_  cr^ 


»  Comparing  the  combined  stress  case  with  the  pure  tension  case  at 
the  same  reliability,  we  obtain  the  fracture  surface 


(W 

lVru, 


,m 

( S0-o  \m 

S0-CT  \ 

+ 

2  u  1  + 

3  U 

cr^. -cr  i 

t  U  / 

t  u, 

m 


(105) 


This  results  in  the  two-dimensional  fracture  diagram  shown  in 
Figure  49  when  S3  is  set  equal  to  <?u>  To  compare  the  approximate 
results  to  the  exact  theory,  we  select  the  parameters  in  both  the¬ 
ories  to  match  the  same  tension  data.  Here,  we  need  only  note 
that  m«k  and  that  the  k's  used  in  Figure  49  correspond  to  those  used 
in  Figure  4& 

We  can  observe  from  Figures  48  49  that  the  approximate  and 

exact  theories  are  identical  for  the  eases  k*»l  and  k*--"  .  Be¬ 
tween  these  values,  the  exact  theory  will  be  found  to  be  the  more 
conservative.  We  shall  determine  the  largest  deviation  between 
the  theories  by  examining  the  case  of  hydrostatic  tension  which 
lies  along  a  45  deg  line. 

In  two  dimensions,  the  approximate  reliability  for  the  hydro¬ 
static  case  is  found  by  setting  Rj  ^ ^ ,  and  »  r^. 
Then,  comparing  this  reliability  with  that  of  the  tension  case, 
we  f ind 

-  (1/2)  l/m  (106) 

t /approx 
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0 


.2  .3 


Figure  49 


Appruj 


Forming  the  ratio  between  this  expression  and  that  described  by 
Equation  (76)  for  the  exact  case,  we  obtain 

(^/Ct)  approx 
^ah/°t^  exact 

This  ratio  is  tabulated  in  Table  VII  for  different  values  of  m. 
The  largest  value  occurs  at  near  m=3  and  shows  a  deviation 
from  the  exact  theory  of  16,9  percent.  This  discrepancy  is  smal¬ 
ler  at  all  other  m's  and  for  all  other  stress  states. 


VtT  P  (nrt-1) 
2  r(tn 4) 


1/m 


m=k 


(107) 


Table  VII 

COMPARISON  OF  EXACT  AND  APPROXIMATE  COMBINED  STRESS  THEORIES 


m 


Two  Dimensions 

(rr  t  rt  ) 

v  h  t/ approx 


=F 


!  (rt'  /"cT  ) 

1  '  n'  f' 


n  t'exact 


Three  Dimensions 

( rr  /  fT  ^ 

v  h  ty approx 

( ry  /  rt  \ 

^  tv  '  t‘ exact 


1 

1.000 

\ 

1.000 

i 

i  2 

1.154  j 

1.291 

1 

\ 

i  3 

1.16935  i 

1.32635 

i 

i 

i 

i  4 

1.163  j 

1.316 

!  5 

1.153 

1.297 

) 

\ 

i  10 

1.110  ; 

1.215 

1  w 

1 

1.000  ! 

1.000 

Setting  S^aS, 

in  Equation  (105), 

we  obtain 

the  approximate 

ratio  of  hydrostatic  strength  to  tensile  strength  for  equal  relia¬ 
bility;  thus, 

i). 

u  ^approx 


(1/3) 


1/m 


(108) 
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When  cru  =  0,  we  can  compare  this  expression  to  the  exact  ratio 
given  in  Equation  (72). 


(s 


t'  exact 


,<2m+-l\ ' 
i(— > 


(109) 


where  we  have  taken  3k=m  so  that  both  theories  describe  the  same 
tension  behavior.  This  ratio  is  tabulated  in  Table  VII;  it  as¬ 
sumes  its  maximum  value  at  m=3.04.  For  this  value  of  m  the 
deviation  from  the  exact  theory  is  32.6  percent;  furthermore,  this 
discrepancy  persists  over  a  wide  range  of  important  values  of  m. 


5.  OBSERVATIONS  AND  REMARKS 


a.  Experimental  Verification 

1.  The  integrity  of  a  basic  volume  is  not  generally  an  in¬ 
trinsic  property  of  the  material,  but  rather,  a  complicated  com¬ 
bination  of  crack  initiations,  propagations,  and  arrestments.  For 
this  reason,  we  must  try  to  characterize  the  behavior  of  a  partic¬ 
ular  basic  volume.  To  do  this,  nominally  identical  volumes  must 
be  used  for  all  the  tests  designed  to  establish  a  fracture  surface. 

2.  The  theory  assumes  that  a  basic  volume  is  subjected  to 

a  homogeneous  stress  state;  consequently,  its  verification  requires 
that  we  meet  this  condition  experimentally.  Furthermore,  the  load¬ 
ing  on  a  basic  volume  must  be  increased  proportionally  throughout 
the  test.  We  note,  for  example,  that  the  case  of  hydrostatic  ten¬ 
sion  produces  no  shear  stress  unless  the  loading  is  incremented 
sequentially . 

3.  If  the  definition  of  failure  is  taken  as  the  inability  of 
a  basic  volume  to  equilibrate  the  applied  loading,  we  must  take 
pains  to  identify  the  fracture  load.  We  must  not,  for  example, 
allow  a  crack  from  outside  the  basic  volume  to  propagate  into  the 
unit  and  cause  its  failure.  Furthermore,  we  must  not  terminate  our 
test  because  a  fracture  has  occurred  outside  of  the  basic  volume. 

When  the  entire  specimen  is  the  basic  volume, this  latter  problem 
does  not  exist;  if  not,  we  face  a  problem  in  the  interpretation  of 
data  which  has  not  presently  been  solved.  With  low  strength  mate¬ 
rials  it  is  sometimes  possible  to  mend  (cement,  glue)  a  fracture 
outside  of  the  gage  length  and  proceed  with  the  test. 

4.  The  construction  of  a  fracture  diagram  requires  that 
strength  values  be  used  with  the  same  reliability  level.  Often, 
however,  one  finds  that  such  diagrams  plot  the  averap  strength 
values.  Unfortunately,  the  average  strengths  from  tests  using  dif* 
ferent  stress  ratios  do  not  generally  represent  the  same  reliability. 
It  is  permissible,  for  instance,  to  use  the  median  stress  for  each 
type  of  loading. 


5.  The  elimination  of  parasitic  stresses  from  brittle  ten¬ 
sion  specimens  is  one  of  the  most  exacting  problems  in  the  area 
of  material  testing  (Reference  17) .  One  can  anticipate  that  the 
problem  will  be  magnified  for  combined  stress  testing. 

b.  Weakest  Link  Materials 

1.  When  failure  at  any  point  in  a  body  necessarily  consti¬ 
tutes  overall  failure  of  the  body,  the  material  is  classified  as 

a  series  or  weakest  link  material.  Now,  if  we  define  overall  fail¬ 
ure  as  failure  at  a  point,  we  create  in  effect  a  series  material. 

In  this  instance,  the  problem  is  to  find  some  way  of  measuring  the 
first  pointwise  failure.  For  example,  if  the  incipient  mobiliza¬ 
tion  of  the  first  crack  in  a  body  could  be  detected,  we  could  use 
this  feature  as  an  overall  failure  criterion  and  thereby  obtain 
a  series  material. 

Because  the  behavior  of  a  series  material  is  understood,  we 
can  separate  the  volume  and  <  ometry  aspects  from  the  behavior 
of  a  test  specimen  and,  thereuy,  make  it  possible  to  identify  an 
intrinsic  property  of  the  material.  Anticipation  of  this  possibil¬ 
ity  prompted  us  to  single  out  the  volume  ratio  ,w/v  in  all  our 
studies . 

2.  Consider  a  component  which  is  constructed  from  a  series 
material.  To  establish  its  reliability  we  can  divide  it  into  n 
imaginary  basic  volumes  and  recognize  that  the  overall  survival 
of  the  component  requires  the  simultaneous  survival  of  each  basic 
volume.  Thus,  the  reliability  of  the  body,  (l-F^) ,  can  be  written 


i  •  f„  -  fr  <i-Ft>  -  n  o*P 

i“l  i“l 


AV .  1 

*~jt“  8 i ^  1  ***2 *** 3^ !  ( 110) 


where  the  reliability  of  the  basic  volume,  1 - F ^ ,  has  been  taken  from 
Equation  (39).  For  an  infinitesimal  basic  volume  we  obtain, 


1  -  F.  -  exp 


zl 
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3.  A  very  large  basic  volume  of  a  series  material  follows  the 
maximum  stress  theory.  To  see  this  we  can  take  the  limit  of  F  in 
Equation  (67)  as  £V -*  00 .  If  this  is  done  we  observe  that  F  »  0  ' 
when  f  »  0  and  that  F  =  1  when  f  ^  0.  But,  f  becomes  differ¬ 
ent  from  zero  only  when  >  cv,  ;  consequently,  if  the  normal 
stress  in  any  direction  becomes  slightly  greater  than  n  ,  fracture 
is  a  certainty.  This  is,  of  course,  a  statement  of  the  maximum 
stress  theory. 


APPENDIX  II 


BEAM  ANALYSIS  COMPUTER  PROGRAM 


In  tnis  appendix,  the  computer  program  wnich  was  employed 
in  tne  analysis  of  the  beam  prototype  structure  shall  be  des¬ 
cribed.  The  program  combines  the  fracture  algorithm  of  Section  II 
with  the  beam  thermal  and  stress  analyses  and  material  property 
curves  of  Section  IV  into  a  single  package.  A  listing  of  the 
program  as  written  in  Fortran  II  for  the  IBM  7094  is  included  in 
this  appendix.  . 

The  program  consists  of  a  main  program  containing  the  thermal, 
stress  and  statistical  analyses  and  of  function  subprograms  for 
describing  the  functions  T(z,t),  E(T) ,  eT(T) ,  k(t),  m(T),  o  (T) 
a^d  oq(T).  The  program  has  been  specialized  to  some  extent  by 
specifying  the  geometry  of  the  beam  and  its  loading.  However,  the 
temperature  distribution  and  the  material  property  versus  tempera- 
ture  curves  are  entered  into  the  program  via  tables  and  hence  the 
effects  of  variations  in  these  curves  may  be  readily  obtained. 

As  an  aid  in  using  the  program,  a  listing  describing  some 
of  the  more  important  variables  entering  into  the  program  is  also 
included  in  this  appendix. 

Listing  of  Significant  Program  Variables: 

depth  in  beam,  measured  from  the  bottom 

size  of  time  interval  between  computations 

ic^  entry  in  table  of  modulus  of  elasticity 
vs.  temperature 

temperature  corresponding  to  E(I) 

function  subprogram  for  computing  the  tempera¬ 
ture  in  the  beam  at  depth  Z  and  time  TIME 

function  subprogram  for  computing  the  modulus 
of  elasticity  at  the  temperature  T 


DEPTH  (I) 

DTIME 

E(I) 

ET(I) 

FBTEHP  (Z ,  TIME) 

FE(T) 


FSO(T) 

FSU(T) 

FTK(T) 

FTS(T) 

FXM(T) 

HIMP 

NE 

NJ 

NPRINT 

NSO 

NSU 

NTK 

NTS 

PFL(I) 

PINT 

POP 

PRATE 

RISK 

RR(I) 


function  subprogram  for  computing  the  Weibull  parameter 
aQ  at  the  temperature  T 

function  subprogram  for  computing  the  Weibull  parameter 
<?u  at  the  temperature  T 

function  subprogram  for  computing  the  thermal  conductivity 
k  at  the  temperature  T 

function  subprogram  for  computing  the  thermal  strain  e^ 
at  the  temperature  T 

function  subprogram  for  computing  the  Weibull  parameter 
m  at  the  temperature  T 

gross  heat  flux  impingent  upon  the  bottom  of  the  beam 
number  of  entries  in  E(T)  input  table 

number  of  entries  with  respect  to  t  in  T(z,t)  input  table 

number  of  time  increments  between  occurrences  of  detailed 
output  printouts 

number  of  entries  in  oo(T)  input  table 

number  of  entries  in  au(T)  input  table 

number  of  entries  in  k(T)  input  table 

number  of  entries  in  e^('f)  input  table 

probability  of  failure  of  the  ic*1  subvolume  up  to  current 
value  of  time 

initial  load  on  beam 

probability  of  failure  of  entire  beam  up  to  current  value 
of  time. 

rate  at  which  load  increases  after  t  *  TIMEP 

risk  of  rupture  of  entire  beam  up  to  current  value  of  time 

risk  of  rupture  in  the  subvoluroe  at  current  value 

of  time 


RRMAX(I)  maximum  value  of  RR(I)  up  to  current  value  of  time 

RRT.(I)  risk  of  rupture  computed  using  the  values  of  stress, 
temperature,  etc.  at  DEPTH  (I)  and  at  current  value 
of  time 

SMAX(I)  maximum  value  of  STR(I)  up  to  current  value  of  time 
STR(I)  total  stress  acting  at  DEPTH(I) 

STRM(I)  mechanical  component  of  STR(I) 

STRT(I)  thermal  component  of  STR(I) 

TIMEM  maximum  value  of  time  for  which  computations  are  to 
be  carried  out 

TIMEP  value  of  time  at  which  load  begins  to  increase 

TK(I)  ic^  entry  in  table  of  thermal  conductivity  vs. 
temperature 

TKT(I)  temperature  corresponding  to  TK(I) 

TRRM(I)  time  at  which  RRMAX(l)  is  achieved 

TS(I)  i1*^  entry  in  table  of  thermal  strain  vs.  temperature 

TST(I)  temperature  corresponding  to  TS(X) 

TT(J)  value  of  time  corresponding  to  WSO(I) 

TTAB(I,J)  temperature  at  depth  ZT(I)  and  time  TT(J) 

WSO(I)  i^  entry  in  table  of  Weibull  cfQ  vs.  temperature 
WSOT(J)  temperature  corresponding  to  WS0(1) 

WSU(I)  entry  in  table  of  Weibull  ou  vs.  temperature 

WSUT(I)  temperature  corresponding  to  WSU(I) 

WXM(I)  i**^  entry  in  table  of  Weibull  m  vs.  temperature 
WXMT(I)  temperature  corresponding  to  WXM(l) 

?T(I)  ic^  value  of  depth  corresponding  to  TTAB(I,J) 
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PROGRAM  LISTING 


CM  UNIAXIAL  STATISTICAL  STRtNGTH  ANALYSIS  PROGRAM 

DIMENSION  iTAB(lOflOO)*ZT(XO)  «TS(10)  f  iST(10)  *E(10)  fETdO)  «TKdO)  *T 
UT(IO)  ,WXM(IO)  •WXMT(lQ)  »WSU{10)  .WSUTUO)  « WSO  (10 )  ,  WSOT  (IQ  t  •  T  T  <  40 ) 

COMMON  TTAtt,2T,TS«TST,E,ETiTlt»TKT«WXM,WXMT,W5U«W3UT,WSU»WSUT»TT»JJ 

1 ♦DTIM£ 

DIMENSION  RR (20) ♦DEPTH (21) » STK ( 2 1 ) ♦ S TKM ( 21 ) ♦ S TR T ( 2 1 ) « SMAA (21) »HRMA 
IX (20) ♦ TRRM ( 20 ) . PR T ( 2 1 ) ♦ PFl « 20 ) 

DIMENSION  2 S TEP (102) *TSTRA<102> , YOUMO ( 102) 1  TEMP (102 ) 

5  WRITE  OUTPUT  TAPE  6*7 

7  FORMAT  (AflHl ALUMINA  H£AM  UNDtR  MECHANICAL  ♦  THERMAL  LOADING  ) 

READ  INPUT  TAPE  5i10»PINT^PRATEiTIMEP 
10  FORMAT  (3F1O.0) 

WRITE  OUTPUT  TAPE  6^ 12^ P INT  ♦PRATE ♦ T IMtP 
,  12  FORMAT  (6HUPINTs,F10.2»3X«6HPRATts»HU#2^ JX»6HTIMEP=«f5t J) 

READ  INPUT  TAPE  5  » 17«  ST  IMEtDUME  ♦  I  !Mt«*  NPHJNT 
17  FORMAT  ( 3F10 .5 » 1 5 ) 

WRITE  OUTPUT  TAPE  6U9*DTIME%T!MtmNPRINT 

19  FORMAT  (7HQnTIMFa»F5*3«3X«6HTlNENStF7*3«3Xt7HNPRlNT8tI9) 

20  READ  INPUT  T APt  5«2H«NJ»NTStNttNTA»NAMtNSU^NS0 
22  FORMAT  (713) 

WRITE  OUTPUT  TAPE  6 ♦ 24 » N J ♦NTS ♦ NE ♦ NTK ♦NAM^NSU tNSC 
24  FORMAT  (4H0NJs«I3»3Xt4HNTS=»J3»3X»3HNts»13«3X»4HNTK8»I3»4X*4HNXMs* 
IF(NJ)40^4U^30 
30  DO  35  J- 1 ♦ N  j 

READ  INPUT  TAPE  5»32« ( T T AB t I»J) ♦ l3l«8) 

32  FORMAT  (OF10.4) 

35  CONTINUE 

4C  READ  INPUT  TAPE  S.42» (ZT I 1) «l3lvN) 

42  FORMAT  (HF1O.0) 

45  IF (NTS)  55t55»A7 
47  DO  52  ! a 1 ♦ N  T  5 

READ  INPUT  TAPE  5 .49 ♦ TST ( I ) ♦ TS t I » 

49  FORMAT  IMVtM  10.6) 

52  CONTINUE 
55  IF(NE)  65^5*57 
57  DO  62  Tsl^NE 

REA!'  Input  TAPE  S.59.ETU)  «E<I> 

59  FORMAT  (Fit*, 10.2) 

62  CONTINUE 
65  IF(NTK)  75.75,67 
67  DO  72  tsl.NTK 

REAL'  INPUT  TAPt  5f69fHT(I),TK(n 
69  FORMAT  (HU.44P10.6) 

72  CONTINUE 

75  DO  9U  fsitNj 

90  TT<l)=FL0AlFU~l)«.5 

WRITE  OUTPUT  TAPE  6.95»  UT  ( I )  «  1  =  1  »fl) 

95  FORMAT  (5hU2T=  ♦F9.4OF10.9J 
WRITE  OUTPUT  TAPE  6^97 
9?  FCRMAt  (6o  TlMt) 

DO  102  J=1 ♦N J 

WRITE  OUTPUT  TAPF  6. 100*TT  ( j)  .  <  T I  At)  <  I  ♦  J >  ♦  I  »1  *«) 
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100  format  ( 1H  ,F5.2,Fa.2,7Fl0i£> 

102  CONTINUE 

WRITE  OUTPUT  TARE  6«l05 
1 05*  FORMAT  (27MQTEMp  <R)  TSTRA I N  (IN/IN)) 

DO  HO  1  =  1  •  NTS 

WRITE  DU  >  Pl'T  TAPE  6«]f7tT$T(l)tTS(i) 

107  FORMAT  (In  ,F 8. ? , 5X . iPE 14, b) 
liO  CONTINUE 

WRITE  OUTPUT  TAPE  bill5 

3  If  FORMAT  (35M0TEMP  (PI  MOD  OF  ELASTICITY  (P5I)) 

no  uo  i  =  l •  Nt 

WRITE  OUTPUT  tape  6i107iET(I)  iE(I) 

120  CONTI nuf 

WRITE  ntilpuT  TAPE  6,125 

125  FORMAT  (J8M0TE.MP  (R1  THERMAL  CONU  (ttTU/F T«HK«R)) 

HO  130  I=1 »NTK 

WPITP  OUTPUT  Tape  6«107,TKt(I)«tM1) 

1.30  CONTINUE 

If (TIMFM-TT (NJ ill  35,1351 132 
332  TlMEMstY (NJ) 

135  TJMf.s-n^  i^e^ST’MF 

KH=H 
i : =1 
J  J  •  2 

CO  137  1  =  101 

ni.pTM<nso.o*«PL«ATp’ ti-u 

SMASH)  =0. 

RRMAX  ( t  J  =0 • 

PFLd'sO- 
J  J7  THRM  t I ) =0  * 

110  Hi  I  =  liNYM 

READ  INPUT  TART  5,  6‘.>,WXMT  (  ()  «WX» U  > 

WRITE  nu  i  PUT  TAPE  6»Ih0iI«W<mTII1 iliWAMti) 

1  .'*0  FORMAT  (hH'iWTMT  ( ,  U  i2H)3,F  lO./diiMMI  , «3X i HHWXM ( ,  U  ,2M>  =»P6.R > 

U1  CONTINUE 

UO  1  A fj  ;  =  i*NSU 

rf.aii  input  tape  nif'JiwbuTi n iwsum 

WRITE  OUTPUT  T.'PK  fe,UJ»IiWS0n  n  »  liwblM  n 
1r3  FORMAT  (M‘»W5U»  ( ,Il»*H>3»F19«*»  JNIK)  oA,<*MWSU< » 1 1 ,2H> -»F ill, 2,5H (PS 
III) 

14S  CONTINUE 

no  i<*a  i  =  iiN50 

RLA'I  ’UP"T  TAP*  Si'i'FiWbaTI  U  iWbOl  u 
WRITE,  nijtput  TAPE  «,d47itiWV»Tm  i  liWVJUl 
J  A  7  FORMAT  HhUWSOT l » l t ,2H) =»FiO*2,3N(H) OXHHWS0  («ll»2H>a,F10,2, SHIPS 
ill) 

1A8  CONTINUE 

150  TIMF=T  tKHtJTTMT 
J  =  J*l 

IF  ( T 1 MF-T  JM'CM )  1  55,155,5 

155  IF(TIMF“TJMEP)  156,156,157 

156  AL=PINT 

00  TO  15* 
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1*7  AL=F  JNT+PRATfcodlMIr-T  IMtP) 

158  8M=1.1«75oaL 
n=o.t> 

E=C.S 

N  Z  “  3  U 1 
PNS=N5 

PNS’-l.G/p^S 
7STtP  M )=-i.n 
! 

HO  160  J  =  .r'Mi 

160  7STFP  (.»)=ZJ>Tm  J-1 )  *PMS1 
iSTfP  ( ►* 4 )  - 1  ,n 

no  1 7n  j=’ 

ZZs-ZiTE1'  i  J)»LK0.5 
ti.-fmtf  mp  ( /?,  t i rir ) 

TEMP  IJ1  ul  J, 

TSTPA  (  !)sf!5Ml| 

170  VOIUKM  l)=F'fe  (T  •  1 

Ha?STkP(»l^)-iiTF»*(r47-U 

PT=0.0 

RMT=m ,n 

pt=h<>»'>  i  t : j f  r a  :  i )  oyotion  < ;>  +tstha  (iut  o^uumoinz)  > 

RMTsHomv  (  f:>THMU)*>YOUMU(l)'>Zi>TLM  1 1 )  *  I  bTH  A  i  NZ )  «YUUMU  ( NZ )  ^ZSY  tP  ( NZ>  ) 

N=n/-1 

no  '.nil  JspiN 

p Tsp r 3 yra  *.n »yuuhoi  j:  «ii 

Ifln  RMT5HMT*H»>n«r.STRA  i.n  »YUUMO<  J»  wZSTfc.PUi  ttlJooH 
SECMUMs:i>.0«M)<‘,,>«#  VS,U 
no  .190 

DEPK-Ur  I»th t  j>  »nm 
n-i>K**  Tm  I  j ; 

r2=F  HTFIIM  ;t|,  t  (f«t  ) 

.ST 9  t  J )  »-►  Tb  (  r 1 1  «rt  (  r;»)  ♦PT/  5  0"^  .  <M)>  t  73ECMOH 

S T n T  t  J i  s-r  I  $  t?P)  “Ft  H?)  *PT  /  (M92.PJ)  ♦iJfc.POHMT/SECMOM 

i<jn  >trm« j'3Tirw(j»-synr u) 

POO  7Quf  t) i  *m**  i n, n ♦  t f’E ; 

TiaPDT^IP 

rvj.iA'i^tTi-iniopAn, 

MK.A  rtN=-1.iMA,">Pr«f  (TO)  /  3»jl»0  » 

HR  V|s<  l  .  /.  /  /I  J.  *>■>•)}’♦<  T  }«©/»;  /  I6')0. 

MIHpsHCM  In^mkAP 
9tSif.  =  U, 

•is  r*si 

00  2<U 
z=.itprr<(  t, 

TEMP  I  1 1  (i.TJMt) 

Ts  TEM!’  1 1 ) 

XH=:V*i^  f  r  > 

SU=PS-Jf  T) 

SO-sKSOlTt 

IF  (Sim  1)->M4.U  n  )  2  m .215  *2 10 


210  SMAxmssrxm 

215  IF  (STR  ( T ' -bU»  220,220,225 

220  H  R  T  ( 1 1  -  0 » 

GO  TO  2 JU 

225  RRT(I)-=0.1»USTRm-SU)/S0J<»*XM 
230  CONTINUE 

HO  260  1=1*20 

!F(R»T(I»-MRTtI*l))  240, 24*', 2-0 
240  RR(D=RRT(1  +  1) 

GO  TO  ?50 
245  RRiPs'JRt  m 

?50  !F(t'4(T»-4MMAx«in  240 * 260 « 255 
255  RHMAX(T)=RWU j 

T R R M  ( I )  =TI**£ 

C=RRU> 

PFt,  U)  =  l.-tX«»F(-C! 

260  RlSkai;T5<«H|nAXiI) 

POF=l.-EXPF {-RISK ) 

WRITE  n'JTpnT  Tapp  6  *  300  *  T I  ME  *  AL  *  8M  *  TbK  AU «  BSTr 
300  FORMAT  (6MnTlMt=,Fb.2,3HMIN,4X.l3HAPPUc.U  LOAU=,F  6  *  1  *  3MLDS*  4X  ,  1 2MB 
1EWJ  nO”fcNT  =  ,K6.1*6H!N-LBS*4X*lUHUMP  cR  AU=  *  F  T  .  1  *  4MH/FT  ♦  4A  *  t>HOSTR  =  * 
2F8.2*3k*PSH 

WRITE.  nUTaiiT  T  APF  6  *  305  iHf.t.T  IN  *  M  IMP  ,*  l  SK  *  HOF 
305  FORMAT  <15*  Nt  l  HEAT  FUJ*  =  *Fo.  2 »  UNO  I  U/H2<»SEC  ♦  Ax  *  16NGK05S  MEAT  FL 
llJX  =  *F6,2*llHHTU/FT2AS£C,4X*5HHlSX.  =  *lPtl4.8,4x»lOMFAIL  PHUb=*OPFiO. 
24) 

<1=1*4  j-n/NpaiNT 

TF(Kl-*2>  150*150*310 
310  <2=F  1 

WRITE  ngipi.T  Tapp  6,j]5 

315  FORMAT  (#>8Hon£MTM  UN)  TtMP  t R )  STRESS  (PSD  MAX  SIR  (PSD 

1  RR  TEMPORARY  ) 

WRITE  OUtP'.'T  TaPF  6*320 

320  FORMAT  tlri  »T0*»50HRR  LAYER  RK  lAYt«  MAX  TIME  (MJN)  Mil  PR 
IOB  ) 

1=0 

330  I=t*l 

IFU-2U  3i5.335.3hO 

335  WHITE  nuTP"T  T»PF  6 . 340 .DEPTH U ) , T tMP ( D « STH ( D , SM«X ( I ) ,KH T ( D 
340  FORMAT  (l.w  ,F7^3*7X»F7.2.6A*Fa.i»M*F8,l«4X«lPEl2.6) 

IFU-20)  300.350,330 

350  WRITE  rilTPuT  Tapr  6*355»«R  t  D  *MR«ax  U  )« )HKH(D  *PfL  ( 1) 

355  FORMAT  DM  ,69X  ,  IPF.lii,6*3X*  lPtl2.6»5**0PP6.3,4X,AF  10,8) 

GO  TO  331' 

360  WRITE  nuTpt'T  TApP  6*365 

365  FORMAT  (54H0PEPTW  (IN)  ME  CM  STKtSS  (PSD  THEWMAL  SIMtSS  (PSD) 
DO  375  1=1*21 

WRITE  pUtP"t  Tapp  fc,J70»CFPTMII)  «STRM(D  »STRT(  D 
370  FORMAT  DM  ,F 5 * 3 *  1GX , Ffl.  1  *  15* »F 8. 1) 

375  CONTINi't 
GO  TO  150 
CNO 


CBT  BEAN  TFMvEmatuhL'  0 f  S 1  P  I BllT  ION  FOnOUON 
FUNCTION  FBTFMh (ZiTIMf) 

OlMtNSTON  lTf.b«i^,lOU),ZT(lC)<TS(10)i1ST{i(j)»t(10)»ET(IU>tTKtlU)»T 
IK  T  C  1C)  .►*►»  ( 1  ^ )  «w*nt  1 10)  ,wSu  ( 1  0)  ,wSin  1 1  0 )  « wSO  ( 1 0)  .  WSO  T  U  0)  ,  I  T(  40 ) 
COMMON  TTAM,ZT»TSiTST,E,f  T  ,  U ,  U  T  » t»AM«  wAMT,WSli«WSUT  « wSIm  wbG  T  »  TT  »JJ 

j  j-  a 

10  IF  (  1  JMF-T  T(  J  .)  J  3COF«2t» 

20  JJ  =  J..'*1 
00  TO  10 
30  1 1 -V 

2*  IF '/-/ T U 1 ' )  hS,uS,40 
AO  11=11*1 
GO  10  35 

it*  TIJ=rTAt  r  I  l,  -j) 

„  TI1J=1TAN(I1-1,JJ) 

TIJ1=TTAM(|I,JJ-11 

TZI  =  TT1.|)*(T*“»--T1  UJ-m«trU.)-Ul.ll)/U.s 
TZ=nji*<Tim  -n  uj-ii  i#n  i  j-i  ijd/o.b 

FUTf,MP?*T/l'*>i?-/r  *n-l)  )»(TF-T21)  /  UT  (  u  )-^T(  II-l  J  )*4i>a.t> 

RETURN 
E  N II 

CFTS  ThEKOAl  STRAIN  A?  A  FUNCTION  UF  I EMF t** 1 OKt 
FUNCTION  FTb'Tl 

OlMFNSTON  ITAhl  ln«  100)  »ZTUO)  «l3»i0;  »’STUU)  <E(iiU)  ttTUOl  *TMiU)  »T 
iKTMfi;  .k.TMlin;  ,uy»Tf  jnl  »U.S0(1D)  ,WSUT  (it))  »W3UtlO)  ♦  WSO’Mlt))  »  T  T  t4U) 
COMMON  TTA^,,,T,T*'1>T*>f  ,L»CT.rF»U  !  «w.»M»WAMT  .wS0,W6U  I  ,wr»W»wSl>T  ,TY  ,  JJ 
1 1 UT IMt 
1  =  1 

10  IFiT-TSTU^l)  )^0,2r,lS 
lb  1=1*1 
GO  TO  10 

20  FYS3Tstn*tr-TSTi:n«tTS(!*i)-isa»>MrsMi*n-Tsr(n) 

RLTUKN 

ENO 


OFF  MO!10t.05  ()F  El.  AST  |  C  t  TY  A  FUNCTION  UF  TtMRtttAUJHt 
FUMCTIrtt,  Ft  1*1 

Dl«F«S ’OR  t  T*'H  !).’•)  00 »  .ZTUOI  •  131  iU»  »T3T  110)  «t  UU)  »FT  UOi  «  IK  UOj  ,  f 
iKTt  ’0)  .w*M<  p.)  .«’»,♦  t  l o ’  »V.V-  { 10)  *wSUT<  iUI  mSO  { 10)  »W3GT  UOt  *  T  T  « t*(;  1 

common  r r am, /r, st ,c,n »tr , if. t *w4H,wAMT.wsu,wbuT,wsu»wsti i .! t > jj 
ItflriME 
!  =  1 

in  ! F  (  T-t T  C  !♦  i )  )  ?o,;os  l*. 

GU  To  10 

,10  FE=r 1 1 ' *•  r-c* i 
BtT'JMN 
EMI 
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CF'fK  THERMAL  COMO'ICTIVITY  AS  A  MJNCTIUN  OP  TEMPtR ATHRE 
FUNCTION  F  T K  t  T ) 

D I  MEN  >  TON  IT'13(  ITilTIO)  -.ZT  (liD  ’  i  SUO)  »  15MIO)  { \  (i)  (l.fl)  1 1IC(  if))  «T 
IK  r  (1U)  .viXMllO)  4  W  X  H  T  do  j  ,WS'-»  (1U)  ,WoUTUQ)  ,  W50  ( 1 0  /  4  WSU  T  U  U  )  1  I  T  ( 40  ) 
COMMON  TTAn,/TiT.S,TST,E,F  T  t  f  K  »  TK  T  ♦  wXM  t  WAMT  «  WS'J  1  wSUT  4  WSU  4  w  Su  T  4  T  T  ,  J  J 
1 4  OT I  ML 
1  =  1 

i.0  I F  <  t-Tk  rtl  +  p  )  ?0,;>0,l5 
15  !  =  I  +  i 
GO  TO  10 

20  FTK  =  TK  (  I'  ♦<  T-T*  T  ( T  M  *  (  TK  (t  +  1)  -  U  <  U)  /  tTKT  U  +  li  -TU  It )  ' 

return 

END 


CFXM  WEI *3‘>L'_  "  «S  A  FUNCTION  OF  TEMPERATURE 
,  FUNCTION  F*M{T> 

Dt  MS  TON  TTAHO.04  ini))  4Zf  (10)  4T.S(1D)  .  TST  f  10)  ,r;<!U)  4f  T  (ill)  4  rx.  (10)  »T 
IKT  ( 111)  ,  4<;m  J)  ,WTMT(10I  4  WSU  tl'))  4WSUT  f  10  5  ,WSl)  1 10)  ,W50T  (  10)  ,TT  UO) 
COMMON  TTA",  ZT«  T$4  fST,C»I  N  r*4TKT  »wXM»WAMT4WSUf  w;,UT»WSU«WoUT»  T  T  iJJ 
1 4  DTI  ME 
1  =  1 

10  !K(T-wXMT (!♦)))  ?0 420,15 
15  1  =  1  +  1 
GO  TO  10 

»:0  FXM-SWX'M  I )  +  (  T-WXMT  (  I )  )«(WXM  ( 1  +  1 )  -W AM  ( 1 )  )  /  ( WXMT  ( I  +  1 )  -WAMT  (  1)  ) 

RETURN 

EN!) 


CFSU  WEIMUU.  sij  AS  a  FUNCTION  OF  TEMPERA  I  URE 
FUNCTION  FbU(T) 

niMENSTON  (  T  > 6  ( 10  4  1.00)  * Z  T  HOI  *  TS<  10!  4)51  (10)  »EU0)  if,!  (ID)  4l  K  (2,0)  ,T 

IKT  (10)  4  W  X  M  ( 1  '• )  ,  W  V  M  T  { 1 0 )  4  WSU  ( 1  J)  4WSU  ft  IQ)  4  WSO  ( 10 )  1 WSU  T  ( 1 0  I  ,  T  T  1 40 ) 
COMMON  TTAh,  7T,TS4TST  ,E,ET*TKi  TK  1 4 W AM, WAMT  4WSU4 WSUV  4 WSU 4 "50 T  , TT  ,JJ 
1 4  DTI  ME 
1  =  1 

\0  JFIT-WSUT  ( l  +  l)  )  90  4  (?0, 15 
IS  1=1+1 

GO  TO  ,  0 

20  FSU='W5<i  ( l )  +  (  T-wS"T  (  1)  )  *>(  WSU  (1  +  D-.XS  J  1 1)  )  '  t/lSUT  U  +  l )  -W5UT  ( 1 )  ) 

RETURN 

EN'J 

CFSO  WEtOotl  SO  AS  «  FUNCTION  OP  UMPlKiUuwt 
FUNCTION  FSO ( T ’ 

dimension  1  tab <  lo,  inu)  ,zr ( mi ,  tsao)  4 1  st  t  iu)  4tuu)  ,f»T « 10)  ,u  ao) ,  r 

IKT ( 10)  ,RX  R 10)  4WXMT ( 10)  ,W3U ( 10)  4WSUT ( 10) 4WSO (10) 4WS0T  t  un  ITT (40) 
COMMON  TTAH4ZT4TS4TST  »E,£T  4TK4  TKT  4WXM4WAHT  ,WSU»  WSUT  ,  'W  SU  4  WSU  T  4  TT  4  J  J 
1 4  DT  tME 
1  =  1 

10  IF  ( T-wSOT  (’•*■!))  ?04?U,I0 
15  1=1+1 
GO  TO  to 

20  FSOswSd  ( I  )  +  (  T-wSMT  ( I ))  »  { WSU  ( l  +  1 )  -RSO  ( 1)  )  t  (w$OT  (M)-«bUltlH 
RETURN 
END 
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APPENDIX  III 

DISK  ANALYSIS  COMPUTER  PROGRAM 


1 


i 


K 


The  computer  program  used  in  the  analysis  of  the  disk  pro¬ 
totype  will  be  described  in  this  appendix.  The  program  is 
similar  to  the  one  employed  in  the  beam  analysis  -  the  greatest 
difference  being  due  to  the  fact  that  the  disk  is  under  a  biaxial 
state  of  stress  while  tne  beam  was  under  only  a  uniaxial  state 
of  stress.  Again,  the  fracture  algorithm  of  Section  II  is  com¬ 
bined  with  the  thermal  and  stress  analyses  of  the  disk  and  the 
material  property  curves  of  Section  IV.  A  listing  of  this  pro¬ 
gram  as  written  in  Fortran  II  for  the  IBM  7094  is  included  in 
this  appendix. 


The  program  is  composed  of  a  main  program  consisting  of 
the  thermal,  stress  and  statistical  analyses  aid  of  function 
subprograms  for  describing  the  functions  T(r*t),  E(T),  8^(1), v(T), 


fc(T),  m(T) ,  ou(T)  and  ag(T).  The  dimensions  of  the  disk  along 


with  an  arbitrary  applied  uniform  stress  at  the  outside  radius 
have  been  incorporated  into  the  program.  The  temperature  distri¬ 
bution  and  the  material  property  versus  temperature  curves  are 
entered  into  the  program  through  the  use  of  tables  and  hence  the 
effects  of  variations  in  these  curves  may  be  readily  obtained. 


As  an  aid  in  using  Che  program,  a  listing  describing  some 
of  the  more  important  variables  entering  into  the  program  has 
bean  included  in  this  appendix. 

Listing  of  Significant  Program  Variables: 


DTIMS 

S(I) 


ET  (I) 


size  of  time  interval  between  computations 

i  entry  in  tabic  of  modulus  of  elasticity  vs. 
temperature 

temperature  corresponding  to  E(I) 


FDTEM3? 

(R, TIME) 


function  subprogram  for  computing  the  temperature  in 


the  disk  of  radius  R  and  time  TIME 


Pressing  page  blank 
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SQf 


iSii 


FE(T)  function  subprogram  for  computing  the  modulus  of 
elasticity  at  the  temperature  T 

FNU(T)  function  subprogram  for  computing  Poisson's  ratio  at 
the  temperature  T 

FSO(T)  function  subprogram  for  computing  the  Weibull  parameter 
oq  at  the  temperature  T 

FSU(T)  function  subprogram  for  computing  the  Weibull  parameter 
at  the  temperature  T 

FTK(T)  function  subprogram  for  computing  the  thermal  conductivity 
k  at  the  temperature  T 

FTS(T)  function  subprogram  for  computing  the  thermal  strain 
at  the  temperature  T 

FXM(T)  function  subprogram  for  computing  the  Weibull  parameters 
ra  at  the  temperature  T 

HIMP  gross  heat  flux  impingent  upon  the  inside  radius  of 
the  disk 

NE  number  of  entries  in  E(T)  input  table 

NJ  number  of  entries  with  respect  to  t  in  T(r,t)  input  table 

NNU  number  of  entries  in  v(T)  input  table 

NPRINT  number  of  time  increments  between  occurrences  of  detailed 

output  printouts 

•MSG  number  of  entries  in  <?0(T)  input  table 

NSU  number  of  entries  in  o'  (T)  input  table 

NTK  number  of  entries  in  k(T)  input  table 

NTS  number  of  entries  in  e^(T)  input  table 

NXM  number  of  entries  in  m(T)  input  table 

PF&(1)  probability  of  failure  of  the  i1^  ring  type  subvolume 
up  to  current  value  of  time 
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POF  probability  of  failure  of  entire  disk  up  to  current 
value  of  time 

RAD  (I)  i1"*1  radius  for  ring  sub  volume 

RISK  risk  of  rupture  of  entire  disk  up  to  current  value  of  time 

RRR(I)  radial  risk  of  rupture  of  i*-^  ring  at  current  value  of  time 

RRRMAX(I)  maximum  value  of  RRR(I)  up  to  current  time 

RRTR(I)  radial  risk  of  rupture  computed  using  stresses,  temp¬ 
eratures,  etc.  at  RAD (I) 

RRTT(l)  circumferential  risk  of  rupture  computed  using  stresses, 
temperatures,  etc.  at  RAD (I) 

RT(1)  i1-*1  radius  corresponding  to  TTAB(I,J) 

S  uniform  tensile  stress  applied  at  outside  radius  of  d?.sk 

SRMAX(I)  maximum  value  of  STRR(I)  up  to  current  time 
STMAX(I)  maximum  value  of  STRT(X)  up  to  current  time 
STRR(I)  total  radial  stress  at  RAD(l) 

STRRM(I)  mechanical  component  of  STRR(I) 

STRRT(I)  thermal  component  of  STRR(I) 

STRT(I)  total  circumferential  stress  at  RAD(I) 

STRTM(I)  mechanical  component  of  STRT(I) 

STRTT(I)  thermal  component  of  STRT(I) 

TXMEM  maximum  value  of  time  for  which  computations  are  to 
be  carried  out 

TK(1)  ic^  entry  in  table  of  thermal  conductivity  vs. 
temperature 

TKT(I)  temperature  corresponding  to  TK(1) 

TRMAX(I)  time  at  which  RRRMAX(I)  occurred 

TS(I)  i  entry  in  table  of  thermal  strain  vs.  temperature 
TST(l)  temperature  corresponding  to  TS(1) 


TT(J)  value  of  time  corresponding  to  TTAB(I,J) 

TTAB ( I, J) temperature  at  radius  RAD(I)  and  time  TT(J) 

TTMAX(I)  time  at  which  RRTMAX(I)  occurred 

WSO(I)  i^  entry  in  table  of  Weibull  oq  vs.  temperature 

WSOT(I)  temperature  corresponding  to  WSO(I) 

WSU(I)  i*”*1  entry  in  table  of  Weibull  vs.  temperature 
WSUT(I)  temperature  corresponding  to  WSU(I) 

WXM  iL  entry  in  table  of  Weibull  m  vs.  temperature 

WXMT(I)  temperature  corresponding  to  WXM(I) 

XNU(I)  entry  in  table  of  Poisson's  ratio  vs.  temperature 

XNUT(I)  temperature  corresponding  to  XNU(X) 
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PROGRAM  LISTING 


CM  B I  AX  I  At  STATISTICAL  STRENGTH  AN AL  75 1 5  PROGRAM 

DIMENSION  TTAB<12tlOO)  «RT  (l?J  ,TS(1U)  »TST  (10)  ,£(10)  »ET(1U)  »H(10)  »T 
1KT (10) ,XNIH10J tXMUT (10) iWXM(id) ,WXM|  110) »WSU(10) «W<>Ut(10) iWSU(IQ) , 
2WS0T(lnj  ,Tf  (40) 

COMMON  TTAn,RTiTS,TST,E,F  TiTttiTKT«XNUiXNUT,WXM,WAMT,WSU,WSUTtWSQ,W 
1S0T,TT,J.J 

DIMENSION  WA  1(70)  »STRR  (70)  '.oIKU  )U)  iSrf(KM(7U)  ,STKRT  (70)  OIRTM(70)  » 
1STRTT  (70)  ,5A'1A»  <  70)  10YMAX i7O),KKiKt7O)»KKTT(7U)*HRH(70)t«KT(/0)» 
2RRRMAX ( 7M) ,RRT*AV (70) ,TRMAX(70) , TTMAX (70) ,  PF  H (70) 

DIMENSION  »STEe (52) ,P0IS5 (52) iTSTHA  <52)  ,  YUUMt)<52)  ,  TEMP  (  7U)  ,RA  (52)  t 
1TPA(5Z,30)  TP A (52,3,3)  iT0IPA(3i3)  , 1NIPA(3«3> 

5  WRITE  0(i ) PUT  Tape  fc.,7 

7  FORMAT  (4,jH],AknMINA  UISK  UNDtR  MtChAN l CAL  +  THERMAL  LOADING  ) 

READ  INPUT  tape  s5,10,S,DTIME»nM£.M,NPKlNT 
10  FORMAT  <3F10.5«I5) 

WRITt  nuTPUT  TAPr  ft , 15 , S , D 1  I  ML , T I *t M i NRK I N T 
15  FORMAT  (  JH«S5tE9.Z,3X,6HDTIMta»Fb.  3«  JX16MT IMEM=,r7,3,3A,  THNPK INT=t 
113) 

20  READ  iMpt'T  TARE  S,22,,nJoNTS,MNU,NE»N)X.,NXM,NSU,NSO 
22  FORMAT  (SIS) 

WRITE  OUTPUT  TAPE  6 »24 »N J , NTS , NNUi Nt »NT* »NXM, NSUi NSO 
24  FORMAT  (4H0NJ=,n,3Xi4HNTS=«  t3»3*^4HNNu=» I3»3X,3MNt=»I3«3Xt4HNTKs» 
H3«3XiitHNXMB.M,4X,4HNSUa,  I3i  JX  ,4rtN5U  =  ,  1 3) 

IF  (NJ)  AO  1411,30 
30  DO  35  Jsl,Nj 

READ  IMpt'T  TAPE  5 »32» < TTaB « I ♦ J ) ♦ i=l » 12 ) 

32  FORMAT  (12E6U) 

35  CONTINUE 

40  READ  IMpt'T  TAPE  5,42i  <RT  < I)  ,1  =  1,12) 

42  FORMAT  (12EC.3) 

45  IF (NTS) 50,50,47 
47  DQ  49  1=1, NTS 

READ  IMP"!  TAPE  5 i4A , T ST <t) ♦ TS  ( I> 

4 «  FORMAT  (f tU.4iFi0.4) 

49  C0NTIK"E 

50  IF (N^yt 55,55,51 

51  00  53  Tsl,NN!i 

RE  All  INP'ij  7  Apt  c,AA,XNUTU)  iXt'UU  ) 

53  CONTINUE 
55  IF (NE ) A5 <45,47 
5T  DO  62  T=1,nF 

READ  INP'.'T  TAPE  5»4R.eT(!)  »EU) 

59  FORMAT  (FK)»4,Fl(l,2) 

62  CONTINUE 

65  IF (NTX) 75, 75.6T 

67  DO  1?.  T  =  UNTX 

READ  Input  TAPE  5 , fe9 , TNT ( 1 1 , T X (  I ) 

69  FORMAT (F 10.4, FIO. 6) 

72  CONTINUE 
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75  00  90  isltNj 
OO  TT (I ) =FUOA  Tp  ( I  —  1 > 

WRITE  niiipiiT  TAP!:  b . 95 « (RT  U  )  . 1  =  1 » 12 ) 

95  FORMAT  <5H0*T  =  »  HF9 . 3) 

WRITE  nu^P'iT  TAPE  fa,97 
07  FORMAT  <5H  r i **e. ) 

00  102  J= 1 . N J 

WRITE  OUTPUT  TAPE  (■ . 1 00 ,  TT  ( J )  .  <  T  T  At  ( 1 . J )  . 1  =  1 .  1? ) 
100  FORMA T  f  lH  .F*»,2.12F9.1) 

10?  CONTINUE 

WRITE  ml.ipi»T  TAPp  f, ,  1 G 5 

109  FORMA T  (  ^ rnuT^MP  (R)  TSTRAIN  I  I N/ 1 N ) ) 

DC  llu  I=1iNTS 

WRITE  n(jipir|  TaPE  6, 1C7*  TS T  ( I )  .  TS>  ( I ) 

107  FORMAT  1 1H  .Fo.2.5x.l?n4*0) 

110  CONTINUE 

WRITE  Mi,  |  p i it  TAPE  6,112 
11?  FORMAT  <?E»Hl'TrMM  fR)  P0IS50N  RATIO) 

DO  115  !=1»NNU 

WRITt  nli  I  ptiT  T  APF  6.1l<*.XNyi  U>  .XNUll) 

114  FORMAT  ( 1H  .rP.*.ftX.F6,4) 

116  CONTIM'E 

WRITE  Ml'T P u T  TAPr  *.119 

119  F0RMA1  (3<9H0TEMP(R)  MOD  QE  tLASTlUTT  <P5I>  > 

DC  120  1=1. NF 

WRITE  "input  TAp F  6*iOMHl)«Hn 

120  CONTINUE 

WRITE  ni.rpul  f Apr  a,U5 

125  FORMAT  f?**HPTFMP  ( R  )  THERMAL  LONU  tb  I  U/M  <>Hk4R  ) ) 
DO  13L  1=1. NT* 

WRITE  OUTPUT  Tapr  ft. 107.TKT ( l >  .T*  1 1  ) 

ur  continue 

IF  (TIMFM-TI  <Nj)  )  U9.i35.13r. 

132  T I M K * s T T  ( N J ) 

135  TI«E=-«'TIME 
j=n 
K2=G 
11  =  1 
JJ=? 

LL  =0 

RAO ( l ) =.9 
DC  137  I=2.7n 

137  flAfl  1 1)  sSQim  tRAlM  1-1  )»«»*♦.  12fr«Ht>9) 

HO  1 39  1=1.70 
ARM AX <  T ) =0. 

3TMAX(T)=0. 

,  RRRMAX  (  I  I  sll, 

RR  TMA  X  ( I )  =U» 

TRMAX  (  n=0. 

TTH'XU)=0. 
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139  PFP ( 1 1 sQ • 

PO  Hi  1  =  1«N*M 

REAL'  InP'it  TAPE  (1)  ,W*I- I  1) 

wKirt  ('iiTP"T  r apf  a.ho,!.wx-m<:) 

J  4O  FORMAT  HHOwX‘"T  ( ,  t  0*t4NWXM(tll»ZH)s»K6«4) 

Hi  roNiiMit 

UU  145)  I=1»NSU 

o F. A l-  iMpnr  TAF't  SH^iWbUI  U  )  iWbu  ‘  l ) 

WRITF  f'ljrpi'l  TM'F  fr«l43,l»W5UT  U  )  «1  iWbUHl 

H5  FoP.r,ATiAH0W5'if(,ji,?H)=,FH'.A<^H(F)<5'X,4hNSU(»lli2h)=»F’lU.2iSH(p51 

'  D) 

14?  CQNTTM"F 

PO  144  ] - i ♦ N  S  P 

RE.AD  lMf  I'T  TApF  M  )  iWSC  «1 ) 

W R I T f  n U  T p  1 1 T  TAPf  4*  l4t> « T » ''"SO  T  1 1 ) » i • W50  •  1) 

144  FORMA  1  I  fHOwSoT  < ,  »  1  ,?H)  =,F1C.2  »  JH  <K)  t^X^4MWS0Ull»iiH)s»^lU.2»5'HtP51 

m 

I4C  CONTTNnf 

15.0  T  I  Ml  =  T  Tl“t  ♦PTIMt 
LL=l  1*1 

IF  M  !KP-tih[m)  J b5> •  1  Sb •  ^ 


lb?  P  =  S«.i'* 

P5TFP (l)sQ.b 

NRSbl 

NT=?n 

R5TKF  UKJsi.n 

PMTsNT 

M-MT 

HO  140  JS»«N 

HO  RSTf.P  (J)  S^srrpu-;  H(?.b/PNTJ 
H=.2b 

no  1  'TO  J  *  1 »  N  T 

TR='RSTF.°  (  J*n  -r.ATFP  t  J))  /?  •  ('.♦MSTtP  «  J) 

Y—FU ThMP  •  Tw  »  T  l Mf. ) 

Tt'MF  «J1  s-1 
VO'.iMCK.t)  -fr  IT) 

TSTRf  U)SF'f  (  T  ) 

P01SS(.n=FNum 

l?c  rcF'Tirtnf 

N  =  f  T 

tJO  1  lb  J  =  1‘M 

Ft  I  J*l)  rl  .11- fRSTfP  I  J>  /  P  5,  r  M  {.'♦))  )  00 1’ 

TP  A  (  Jf  \  ,  1 ,  n  =  (*•  O’  F  ('  I  J*  1  )  /  kM  I*'  •  J  )  >  t  i  •!'"  I  i  .U*PUI  5S  I  J)  )/2«U0Kaij*>) 

TP  A  I  J*  I  «  i  »2)  -K-Sl  FF  (.1*1 )  ®  ( 1  •  (i-K'Iba  1  J)  ''‘'i'l  ( J*U  /  I  YOUMO  I  J)  oh<»2,U) 

rp  A  I  ,«♦!  ,  1 ,3)  =  (  1  .*♦(*'' 1  5S  (J  )  I  /^.IIOTSTH*  l  J)  “MSTtP  J  J*ii  »PAi  J*D 
TFAIJ*)  «/«)):(  vO'.'T’i;  (  J  )  <>KrPt  1  1 )  I  /  l  c  .  «*R  S>  f  tP  I J  )  ) 

TP  A  I  J*1  , p.<*}  =1  .(j~  M  .  (1-POT  St  I  J  >  )<>WA'J*))/<’*0 
TPA ( j*1  *2  1 .1)  (  J)  *’l  *  T  K*i  l  j  I  °KA  l  )  ✓<  »C 

TPA«.l*l,  i,U=U*0 
TPM  J*1  ,  W)  sO.o 


17?  TPA<J*l,i,3)=1.0 
CV1PM1  ,1,11=1.0 
OVIpAtl. 1,21=0,0 
QVIPAd  ,],.•»)  =0.0 
OVIPAtl  ,2,l)s0»0 
OVIPA (1*2,2) =1*0 
OVIPAtl,  2,i)=0.0 
OVIPA  0 ,3,1 >=o.o 
OVIPAtl, 3, if!=0.0 
OVIPAtl, 3, 31=1.0 
CO  1*0  t=l»3 
DO  1*0  1=1.3 

160  OV  tPA  t2,R  ,U)  =TPA  (t,lf  U  5 
M=NT+1 

IPtN-31  1'?*,  161,1*? 

16?  DO  190  J  =  J,D 
DO  190  L=l,3 
UO  190  f'i  =  l ,  3 

190  OVJPA  t.l«l«»*)*TPA  t  J.l,I)*CVlPAt  J-1»1i*»*TPa(J*1,Z>  *OVIPA  l  J*1  »i?  ,*)  ♦ 
11PAU,L,35«*0VI**A(J-1, 3,P) 

Tps  ( P»f>vIP  A  ( 9T*1  ,i*,  3>  1  /OV  I*A  INI  ♦  !  1  ) 

DO  24C  K=l,70 
J=1 

P  1  =P  A  (I  (  ) 

?on  Rtak&iFr « ji 

IF  (Rp-ri)  pr,s.<?io,?ps 

20?  J«J*1 

GO  TO  •»0O 
2i0  CC  m  !'h? 
no  ?  j ?  1*1.3 


21?  TOlPAtl*t»sOVIBAtJ»r,L* 

IF  t  J»1 1  >3'»..?3?,22r 

220  jsj-1 

GO  in  234 
22?  J-J-i 

c F  1 1*. i  .f>«Ri?r*P<j>/4  1]u><'2 

TNIPAM  %l)i=tol/Ke*'rt*<Jl  ><jC!»t-"“<1.0*F>O*SStJl  1/2«0<»RP  11 
YM  P  A  1 1 ,  J* )  =K  i  «>  *  1 » ii“{,n  1 SS  t  J 1  w<*2  •  °PF  )  /  *vOUMUtJ)*>H«>2,Ul 
TKIPAtl  ,  ’<;*(!  .0*'1C!SMJ)  1  /2 . 0°  T  3  T  m  a  t  J )  «K  1 

YfiD'MM  )*  ivQi'N'M  J)  <‘M«KF  1 )  /  f  2 . 0<>H ? I fcK  t  J )  ) 

TNTPM*./)  =5  »tl-«  '  .D-POlSSt  U  )«MF  V/2.0 
Tf'lPf.  (?,  H)  3»^«Yynnnt  J)  wYSTKA  » J )  «*MF  1/2*0 


TfilPA  ( ■»,  i ;  =6,0 

tmita  m,*»=o,0 
TN1PA<3, JJ  =1.0 

tie  230  1  =  10 


no  230  mo 

230  TOIPAH  ,l}«TNI<*Aa.  t)»OVlPA\J«l«I)*lNlPA(L«2)*0VtPMJta.l>* 

lTN>IPA(L,31o0VlPAfj,?,n 

235  0lSPatf»T0tpau,i)*Tt)iPAa»3l 
T  =  FnUMP<R) ,YI«Kl 

5  TPH  l  K  }  s  tTP'*T"Ir‘t*\2,l)«'rO)PAt2,3n/l' 

_ STPT(KI=F£ IT) »DI?P/01»FNUt ‘ ) °S iMk  U ) -F £ <1) "f YS ( T > 
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TP  Os-OV  I  PA  (NT*1, 20)  /0VIPA(MT>iOtl) 
ni5PasTptH»tOTpn(l,jl+TUTPAil»3) 

STRRUKJ  =UPO*TO!Pt  (H,l}  +  Tt;iPA(2,j)  )  /  M 

STRTT  (»  )=FE  <T)  *I)ISPn/f{l4Fiynj(  i  )*5ikrt  (M-PtM  J 

STRPMOt  )sSiRR  (H)-STRRT(K) 

STRTM(K)=51RT(H)-STRTTU5 
?A0  CONTINUE 
250  TO=F  UTfmp  !li,5,  U-'E) 

Tl=FIJTFHP(U,s;>*YtMF) 

TGRAL’s'Tl-'U)  PXAO. 

HE  A I  I  Ns-l  GH  Atl#F  TP<T0)/3600» 

HP AD=  ( \  .  I  \  f/10«<»«9]  <>(T0<><>4)  /3b DO* 

himp=nfatin+hhau 
P I SK  =0  * 

RISKRSO, 

RISK  1=0. 

00  30C  1=1*70 
TR=RAU ( I ) 

TEMP  {  n»MJlt:«P  (TO, TIME) 

TaTtMP ( I ) 

XM=f  XNtT) 

SOSFSO(T) 

SO=FSO(T) 

IF  <STRB  t  n-SOMAX  rn  Uf.5«‘fc5«260 

?6G  SPMAX ( ! )=STRR 4 I) 

?bS  I F  ( S  TRW  1 1  )-SU)*7n,;??0«ii75 
?7P  PR  TP  rpfC, 

GO  TO 

?7b  pptp i  (strp i n-su)  /soh>»xm 

270  IF  (STKTm-STMAXU!  )  t6$*<-85*ciJ0 
?6C  STMAXITl&STPTtt) 

?-85  IF  (SIM 
2*10  PRTT  ( 1 )  asti. 

GO  TO  3GU 

?9S  RHTT  1 1 )  *.090597$*$*  C ($TRT  ( I )  -6U)  /5U)**XM 
300  COST  1 NoF 

DO  350  isi.o1) 

IF  (PPT>M  \)  ->fu;h  ( ) )  3lC’210* 318 
3 1 C  RRR(l)=P»Twtr*l) 

no  to  n<) 

31$  nRRU  JrpPTw(I) 

32C  IF  (RPP  1 1  >  -KpuHAy  r  J )  )  3?fl*  ■»?♦**  325 
32$  RRPMAX ( ] i =mrr  1 1 ) 

TftMAX ( T )at Ihf 
32e  CR=PH»maM!5 

IF  <RPir  ( i > -fptt  i  1*1 )) 33C , 3i$ *  3 3$ 

330  PftTtl)=PWTUW) 

GO  TO  3MI 
335  Pfc  T  ( I )  aRR  T  I  vi) 

3A0  IF(HPH  D-RPTMA/  (!)  )  3<,e.3<»P,3<*$ 

3A5  RHTMAXf T)=RRt(IJ 

ttmax(»)=timf: 
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246  CT=RKTMU  ( I ) 

PFR (I ) s] .-tXPF  t-CP-CT) 

RISKK  =  hISKk+«RRM£iX  {  I ) 

350  RISICT=P!SKIfPRJMAXn) 

R I 5<SP ! SKR*R  t S* T 
POFR=l.-EXRFt-Mir>KRl 
POFT=l,-EX^F  (“RI'XKT) 

FCF  =  1  * -t* l-RISk"  ) 

WRIU  niiipuT  TaPp  6».3r,U,T]ME*:>,Tu»Tt,kMj 
360  FOPMAT  ( GHO  T  I  TE=U5.2,3HM:n,4X»150APPL  1  LU  STKtSS=« F  / .  i » JKM?  I  »4X , 10H 
1 S v J R F  TFMP  =  ,fA.l,lHp,sy.,l(nntMP  GMAt)=,P7*lJ 
WPiU  nrrpnT  tapf  e  » 3  f  b  •  h  f  *■•  t  i  n  ,  h  w  a  l>  ,  h  i  m  p 

366  FORMAT  (  NET  hr  AT  F  LU'  =  %  F  6  .  <d , 1 1  nM  U/ H  id^StC  ,  4X  ,  9hhE  A  T  mAU=,P4,£, 

T2*»5FC*Ay  ,if,Mr,RO?S  HEAT  U  ux  =  i  F  fa.  i , 1  lHbTli/c  T2«»&fcl. ) 

WRIU  OiMPu’t  TApt  6,J70,POFR»POF  TiKlSMPUF 
370  FORMAT  USH  R  A  L)  FUL  PpQB=  U  10  . * , 4X ,  U>M  taN  NIL  PK08  =  U  1U .  B  t  4X  »  13H 
iOVIRALl  RISKS,  1PF14.8.19H  OVERALL  FAJU  PKUB^OPHO.e) 

Kl  =  J  +  U.L-\)  I  ,;T 
IF  UMt)  IbC , lbC  *  3'F 
,376  K2=K  1 

WR  I  TP  r'l'TP'iT  T  aP t  6,380 

380  FORMAT  RAHTHS  TEmP  STRUM  STMTM  STRhT  5TR 

ITT  STWK  STRT  SK4AX  SIMA*  ) 

00  3U0  1=1 «7n 

WRITE  '■'MPtiT  TAPF  h  ,  ilpb ,  R  Al1 1 \)  ,  TEMP  1 1  >  ,  STKKM  U  )  ,  STM  T  *  U  J  ,  STHMT  ( 1  >  , 
1STRT  l  (T>  ,6  l  PR  1 1 )  ,8TPT  (I )  ,  SPMAX  *  1 1  »SfMAXU) 

385  FORMA!  t  in  ,FH,4,9Hn,0) 

390  CCNTIM'C 

WRITE  OUTPUT  T APf  6,395 
396  FORMat(?'»mo  pah  «r*r 

WRITE  OUTPUT  TAPF  6,400 
400  FORMAT ( in  «34X*Y6HRRH 

1  PP1MAX  TTMKX  PER  ) 

1=0 

Qf‘FR=U, 

PF«<70;=-1. 

410  1=1*1 

JFIPPkC)  *QPFH»4?f',410,4J;0 
4?0  WRITE  ri.TRUT  TppF  6*425  ,pAl'  <  I  )  »RM  *R  U  )  »MMT  T  ( () 

4?.F  FORMAT  ( 1*  «t  A ,  3 ,  IFF  I  4 . 4, F 14 .6 ) 

IF  U-6<5)43"«43»>0  5" 

430  »F <Ft*K i m 43P,431.m3E 
4J 1  OPTRSO. 

GO  TO  MM 
432  OPFRsPFRlp 


WRITE  %lPiiT  TAPF  8,435»RRR  II  )  «RK  I  U  )  «RHKMAX  (  I) -tRMAX  (  1 )  ,KKT«AXU) 

MTFamn  ,Mr»m 

435  FORMAT  ( 1*’  •35x«I*tl*.ft,2U‘*.*«OPE  f '  *  -  ARt  14.  b  ,0pM,  S ,  lpt  14.6 ) 

GO  TO  4»  1 U 
FNIJ 


MM  T  T  ) 

MKT  RRRMAX  THMAX 
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\v 


u  *  4  * 


CDT  DISK  T^'N^EH ATUHE  DISTRIBUTION  f-UNCTlUN 
FUNCTION  FlJTFMH  (R,tIMn 

Ul^flSTCN  !  TAB  1 1?»  i  no  i  ,RT  (J  ?!  *TS  t  IP)  iT?t  ( 101  *E  (  HU  ,fcT  '  j  0 1  *  U(  10)  »T 
IK T  ( 10)  ,  XMUv  ) r')  ,Ak'UT  ( 10)  tWXN  ( 10 )  ,WXM  UO)  »WbUl  10)  *W5UT  (10)  ,v»5U(  10) 
2WS0TUD}  ,Ti  (40) 

COMMON  TTAM,RT,T5«TST,E»FT«TK tTKTO.NU«XNUT,WXM,WXMr ,wSU,WSUT«WSOiW 
1S0T,TT,jJ 

1C  IF tTlNF-ITt j j  J ) 30*30*20 
20 

GO  TO  Ui 
30  TI-2 

3b  IF (R-RT ( I  1 1 ) 4b«HK«4n 
40  Ilan*j 
GO  TO 

4b  ti j  =  tTf e  !  I  I  ,  JJ) 

TI1JSTTAMU  1-pjj) 

Tl  J1  =  Tt^  ( i  ]  ,jj-]  ) 

TIUl  =  T7ABt  1 I-PJJ-1) 

TRI-TP  Jl  +  (TlNt-TT{Jj-lJ  )<M  Til  J-TllJU 
TR  =  Tl  Jl*  (  T  IMF-1  T<  JJ-1))  MT1  J-TI J1  > 

FOTEMP  =  twi^(s-rt<  1 1-1)  )<MTH-mi>  / (K T  < 1 1 ) -XT U 1-1 ) )  *4b'».6 

RETURN 

ENU 


CFTS  THERMAL  bTMATN  Ab  A  KiNCl  ION  Oh  RN'Rt>ATUHt 
function  F T S ( T > 

OlMtNSTCN  rT&8(l?')OC)«RT(iS),i$tl(J)'1SmO)«C<lU)«tT(lUj«TIUlQ)«T 
IK  T  ( 10)  ,XN'J(10»  ♦X*iUT(10)  «WX*(1U>  »WXMT  UO>  «WSU(10)  i"MJT(iu)  «wSl>  (10>  , 
2WS0T(ln) .Ti ( 4y ) 

COMMON  m*,«T.Yb»1ST,E%F T»T*»HI  i*NU»ANUT»WAM,WAM1  ,W5U»wt>UT  *WS0«W 
ISOTiTT,JJ 

i*r 

10  IF  C  T -  T  SM  1*1)  )  P0,20«  lb 
lb  1=1*1 

GO  Tfi  10 

20  =T?arS»  l  J*iT-?j»l  «  M  )  *MTS<  J*1)-TSU)  )  /  t  T  5  T  INU-TST  U)  ) 

PE  TURN 

F.vo 


CFF  MC»)I>U,<  OF  El  A^T’CITV  Ab  A  F uNC T  l Um  0»>  1  tMRtRA T OHt 
FUVC  rjor.  Ft*  r  *) 

t'lME*'b’r«^  i  i  Ab«  1?,  )o0)  I  (1?)  ,  Ibtlo)  PSP  Ml)  »i(  lu)  UU»  »  U  (10)  *T 
IK  T  ()(*}  OMII  J'’)  •X*:i.*T  ( IC  )  ’XrilOI  i  WXMT  (10)  »WSU(10)  fWSWTUOl  iHbU  UO)  n 
2WS07  tU)  «TT  i  i*U) 

COMMON  TtAH,RT,T<;,TST»EfFT«Tlt  T  »  XNU«  XNUT  *WAH,WXM  I  »  WbU  ,  WbU  T  «  WSO  «  W 
isor,TT,jj 
1  =  1 

iC  IF(T-ETU»m  i'O.i'CPb 
15  1=1*1 
GO  TO  ip 

2C  FE  =  t ( l ) * i T-fc. T (1 ) ) « » E ( !♦ 1) -L ( I ) ) / (t t ( i*j ) -tT ( 1 1) 

RETURN 

ENti 
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CF>ij  P01S5UM  HAt  in  As  A  FUNCTION  OF  UMPEkAt  u*E 
FUNCT Inr.j  FMl(T) 

DlN‘*EftSTON  I  TAb  t  IT  ’  10'))  %RT  ( IC)  »  TS  ( It!)  i  TST  ( 10)  iE  ( 1U)  iF.T  ( 1U!  »  U  m  o)  ,  T 
1KT  ( 10)  ,xvij'  I'D  .X'  UT  (10)  ,WX"UvJ>  ,w^T(10)  ,WbU(lC)  ,wbUT  1 1 0 )  ,  w^O  !  10)  , 
2VS0T(lp) *Tl (40) 

COMMON  TTAM,BT,TS«TSTtE,t.T«T<»llCI  i  XNU ,  XNUT  « WXM,  WXM  r ,  wSU ,  wSUT  ,  *'5>0  » W 
1S0T«TT«JJ 
is) 

io  if  n-xMur  ( i*d  »in,?o,is 
}.?•  i=t*i 

GO  TO  U 

20  FMIsXM*  M  )  *  (  t-xkuT  n  )  )  »(XnM  1*1 )  -A\U(1)  )  /  UMIT  (  l  +  l)  -  XNWTI  i  )  ) 

RETURN 

EN!) 

CFTK  THERMAL  CQNOmC  I  1\'1T  v  AS  A  HjNUlUN  OF  T  L  HP  tK  A  T  )ik  E 
FUNCTION  F l n  t  T  > 

DiMFNST?N  tTau(12»10U)*P?M2)tTMin).TST(in>,t().a)»tTUt)>«T<(10)*T 
IK  T  ()  0)  «  xMJ  t  ID)  «X‘.UT  (  10)  ,V.XM(1D»  ,wx**.l(  10)  iWbUUG)  »W5UT  UU)  ««.'SUUO), 
2WS0TI 1 o) , T I (40) 

C0M*'UN  TT  AFi,VT«TSil  *  »  »E»C  T«TX  »  iK)  *  AM,1 i  x^U T « WXM,  W AM  I ,  W'JjUi wiUT  » ViSOi  W 
1S0T«TT,J.i 
1  =  1 

10  IF  (T-TVT'I*1)  )  2C*«;')»lb 
15  1=1*1 
r,c.  TV  K 

20  ftk  =  tk  ( 1 1 *(t-tk Tun«(TK  u*u-ik  tn  >/(tx tii*d-tk  t  in ) 

RETURN 

EHD 

CFKM  WtlBi'Ll  *  *5  A  FUNCTION  Of  TMVtHMUkt 
function  F*M(T» 

D I  Mf  *.5  T  C*'*  I  T  AO  u  7 « 1  ni!  j  %  T  S «  U* )  •  <  S  T  (101  •  E  f  1  0 )  «  E  T  t  1 0  >  «  U  »  It')  » T 

1KT(  1  Q)  ,  X‘U'H  l'>)  »  XNUT  •  10)  »WA-M  1U)  ,W*M  HO)  ,*bU  ( 10)  *  wbUT  1 10)  i  wbU  ( 1U>  i 

awsnmo)  «tt  (4<j) 

C0MM0N  TTAr»RT»T<IiTS!,EiCT*Ti'»TKT»*Nu«XNUT,VJAM,w»Mt,if,l>U»i'bUT«*J!iiO»W 

1S0T  »TT , jJ 
T  =  1 

TO  IF  (T-wXMt  (  t*T  3  )  ?t,.2'),le, 

A*  1=1*1 
60  U  it 

2C  FXMawX»(i)*(T-WX’«T(n  )<Mv!XMU*l)—FH«i))/(WXHT«Ul)-WXMMU) 

RETURN 

EM) 


cfsu  wti full  su  as  function  of  t t^p tKA tune 

FUNCTION  FSU(Ti 

niMesiST0^  TTAb(l?'inO)  ,RT  ( J  2)  'TS(IU)  « 1ST  (10)  'tl  UM  «t.T  uu>  t  J<  ( 10)  ,  T 

1  ^  T  UOi  »  >ru  1 10)  tX*:UT(  10)  Il7f,']0)  .WX*  I  ( iy)  »*5U  1 10)  »W5UT  < 1UJ  .wbU  ( 10)  « 
2WSOT  l  If  ) ♦ T )  (40) 

COMMON  TTAH,PT*T*.,’rsi  tE,FT*TR*  TIC  1  » >NU« XNUT , WXM . WXM T , WSU»  w$uT  *  wSO «  W 
1 S  U  T  ill  » j  J 
1  =  1 

if  tF!T-^m  tr  +  n  )  rc,?o,i5 

15  1=1+1 

GO  to  -10 

2C  FSiiswSi'  ( I )  +  C  T-WSl-T  ( i )  )  »  ( WS1‘  1 1  +1 )  -wSU  ( i  )  )  /  t  WSliT  ( 1+1)  -wSlJT  i  n  ) 

RETURN 

ENO 


CF 50  WEIbUH  So  AS  a  FUNCTION  OF  TtMRfrKttTUKE 
FUNCTION  FSO(T) 

DlME'‘iT0M  1  Tflb<l<ilr,U)  » H  T  ( 1 2 )  «  T  S  1 1 0 )  «  T  5  T  1 1 0  >  « F.  ( 1 U 1  t  E  T  1 1 0 )  t  T  X  1 1 0 )  .  T 
l  X  T  (10)  ,X»HJ  1 10)  .XNUT  (10)  iVX^  ( )  U)  ,wXMl  (  iu)  ,v»SU  C  J  C)  .WSUT  ( 1U)  ,  wSU  (JO). 
2WS0T 1 10) ♦! I (40) 

COMMON  TTAk  oy.T<;,  TST  ,f  ,fT<TX  ,  K  T  ,  *  MU « X^UT ,  W*M  ,  WXM  I  ,  WSU » wSU  T « WbQ  ♦  W 
1 SUT  ,  TT  ♦  J J 
1  =  1 

10  I F  t  T-wSOT  U  +  l)  1 
15  1=1+1 
GO  TO  ui 

£0  FSO=wSn  (I) + 1  T.v$*T ( !) )  » (USO  < )  +1 )  -w!jU  a  )  )  '  t  WSOT  ( 1  ♦  1 )  -w5U  T  (  n  ) 

RETURN 

FMi 
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